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ABSTRACT 


A  system  that  incorporates  distributed  digital  subarrays  working  cooperatively  as  a  single 
array  can  potentially  increase  the  output  signal-to-noise  ratio  and  provide  better  spatial 
resolution  compared  with  using  the  subarrays  individually.  However,  collectively 
combining  periodic  widely  separated  subarrays  results  in  unacceptable  grating  lobes,  and 
these  lobes  cannot  be  suppressed  using  traditional  windowing  methods.  In  this  research, 
we  focus  on  distributed  subarray  antennas  that  are  comprised  of  subarrays  that  can 
operate  individually  or  collectively.  We  develop  techniques  for  grating  lobe  suppression 
on  both  the  transmitting  and  receiving  sides  of  the  distributed  array  system.  Traditional 
solutions  and  new  methods  are  examined  in  detail  via  numerical  simulation  to  quantify 
the  performance  limitations  when  applied  in  combination.  One  contribution  of  this 
research  is  a  hybrid  approach  that  uses  a  combination  of  suppression  techniques  on  both 
the  transmitting  and  receiving  sides.  Another  contribution  is  the  development  of  new 
receiving  processing  methods  to  suppress  grating  lobes  and  improve  the  signal-to-clutter 
ratio  and  signal-to-interference  ratio.  A  final  contribution  shows  the  relationship  between 
thermal  noise,  array  errors,  and  the  grating  lobe  suppression  effectiveness.  The 
consideration  of  array  errors  addresses  the  issue  of  array  calibration  and  synchronization, 
which  are  critical  concerns  when  multiple  arrays  operate  coherently. 
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EXECUTIVE  SUMMARY 


Complete  digital  control  of  amplitude  and  phase  at  the  element  level  of  an  array  allows 
great  flexibility  in  beamforming.  Modern  radar  and  communications  systems  incorporate 
phased  arrays  with  wider  bandwidths,  allowing  for  the  possibility  that  several  systems  on 
the  same  platform  can  share  arrays.  A  system  that  incorporates  distributed  digital 
subarrays  (DDSAs)  working  cooperatively  as  a  single  array  (thus  forming  an  array  of 
subarrays)  can  potentially  increase  the  output  signal-to-noise  ratio  (SNR)  and  provide 
better  spatial  resolution  compared  with  using  the  subarrays  individually. 

Another  factor  impacting  architecture  is  the  platform  design  philosophy  for 
military  applications,  which  has  changed  dramatically  with  the  advent  of  stealth 
technology  and  requires  reduced  platform  signatures.  Due  to  the  stealth  requirement,  it  is 
difficult  to  find  an  available  area  sufficient  for  a  large  array  on  board  a  ship,  so  it  might 
be  necessary  to  use  several  relatively  small  noncontiguous  (separated)  areas  (subarrays) 
and  then  process  the  received  signal  coherently. 

Traditional  periodically  distributed  subarrays  (subarrays  whose  centers  are 
equally  spaced)  form  a  long  baseline  and  are  capable  of  very  accurate  angular  location  of 
targets.  However,  collectively  combining  periodic  widely  separated  subarrays  results  in 
unacceptable  grating  lobes,  and  these  lobes  cannot  be  suppressed  using  traditional 
windowing  methods.  Grating  lobes  appear  in  the  visible  region  if  the  subarrays  are  nearly 
periodic,  and  they  are  unwanted  because  of  the  ambiguities  that  accompany  them.  Even  if 
the  individual  array  patterns  have  no  grating  lobes,  conventional  beamforming  with 
periodic  DDSAs  will  have  an  output  response  with  grating  lobes,  which  is  unacceptable 
for  most  applications. 

In  this  research,  we  focus  on  distributed  subarray  antennas  that  are  comprised  of 
subarrays  that  can  operate  individually  or  collectively.  It  is  assumed  that  no  grating  lobes 
appear  in  the  visible  region  for  each  subarray  when  scanned.  We  develop  techniques  for 
grating  lobe  suppression  on  both  the  transmitting  and  receiving  sides  of  the  distributed 
array  system  (DAS).  Traditional  solutions  and  new  methods  will  be  examined  in  detail 
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via  numerical  simulation  to  quantify  the  performance  limitations  when  applied  in 
combination. 

In  Figure  1,  we  summarize  the  effectiveness  of  the  conventional  methods  and 
their  combinations  on  grating  lobe  suppression  for  this  specific  arrangement  without 
amplitude  tapering.  The  improvements  shown  are  computed  relative  to  the  periodic 
DDSA.  For  concise  presentation,  we  define  the  shorthand  notation  of  each  method  as 
below: 

P:  periodic  DDSA.  RD:  random  subarray  displacement. 

SR:  sequential  subarray  rotation.  RS:  random  subarray  sizes. 

In  Figure  1,  it  can  be  seen  that  the  combination  of  methods  has  the  greatest 
improvement  relative  to  using  the  conventional  methods  individually.  Flowever,  it  is 
worth  noting  that  an  optimization  process  is  needed  for  each  new  DDSA  arrangement  to 
achieve  maximum  grating  lobe  reduction. 


ca 

"O 


10 

5 

0 

-5 

-10 

-15 

-20 


irrif 


P 

P+RD 

P+SR 

P+RD 

+SR 

RS 

RS+SR 

RS+RD 

RS+RD 

+SR 

■  Largest  GL  (dB) 

-11.03 

-12.33 

-12.72 

-13.95 

-15.96 

-16.22 

-17.95 

-18.95 

V  Improvement  (dB) 

0 

1.3 

1.69 

2.92 

4.93 

5.19 

6.92 

7.92 

Figure  1 .  Summary  of  the  effectiveness  on  grating  lobe  suppression  using 
conventional  approaches  individually  and  in  combination. 


One  contribution  of  this  research  is  a  hybrid  approach  that  uses  a  combination  of 
suppression  techniques  on  both  the  transmitting  and  receiving  sides.  The  result  is  an 
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improved  two-way  pattern  performance.  For  comparison  purposes,  a  two-way  pattern  of 
the  periodic  DDSA  is  shown  in  Figure  2.  The  two-way  pattern  shown  in  Figure  3  is 
generated  by  multiplying  the  transmitting  pattern  (random  subarray  sizes)  by  the 
receiving  pattern  (random  subarray  sizes  and  virtual  filling).  The  side  lobe  level  has  gone 
down  to  less  than  -50  dB  without  affecting  the  mainbeam. 


o  □ 

Figure  2.  Periodic  DDSA  two-way  pattern  of  signals  from  10  and  -25  related  to 

the  broadside. 


Figure  3.  Random  subarray  sizes  DDSA  two-way  pattern  of  signals  from  10  and 

o 

-25  related  to  the  broadside. 
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Another  contribution  is  the  development  of  new  receiving  processing  methods  to 
suppress  grating  lobes  and  improve  the  signal-to-clutter  ratio  (SCR)  and  signal-to- 
interference  ratio  (SIR).  We  propose  virtual  filling  of  the  gaps  between  the  subarrays  to 
eliminate  the  grating  lobes  on  the  receiving  side  so  that  the  response  of  a  single  large 
contiguous  array  is  synthesized.  Therefore  no  grating  lobes  will  appear  as  long  as  element 
spacing  within  all  subarrays  is  less  than  one  half  of  the  wavelength.  Furthermore, 
amplitude  tapering  can  be  applied  to  the  synthesized  array  to  reduce  interference  and 
clutter.  Consider  a  five-subarray  DDSA  with  thirty  elements  in  each  subarray  and  an 
element  spacing  of  0.42/1.  The  subarray  length  is  12.6 A  .  The  gaps  are  also  (arbitrarily) 
set  to  12.6 A.  One  unit  amplitude  signal  is  incident  from  0  with  a  phase  of  7r/5.  A 
second  interference  signal  is  coming  in  at  2.3  with  a  phase  -4tt/5.  In  Figure  4,  the 
average  synthesized  array  response  with  an  SNR  per  element  of  6  dB  (single  snapshot)  in 
the  direction  of  both  signals  is  the  same  as  that  of  a  contiguous  array  as  the  weights  are 
changed  to  scan  the  main  beam  in  a  region  of  direction  cosine  space  (sin#) .  A  20  dB 
Taylor  amplitude  distribution  is  applied.  As  can  be  seen,  the  high  response  of  the 
interfering  signal  that  occurs  at  grating  lobe  locations  has  been  eliminated. 


Figure  4.  Comparison  of  original  and  synthesized  antenna  response  after  virtual 
filling  for  an  element  level  SNR  of  6  dB.  Taylor  amplitude  taper  ( n  =5,  SLL=— 20  dB) 

has  been  applied. 
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A  final  contribution  shows  the  relationship  (and  hence  tradeoffs)  between  thermal 
noise,  array  errors,  and  the  grating  lobe  suppression  effectiveness.  The  consideration  of 
array  errors  addresses  the  issue  of  array  calibration  and  synchronization,  which  are 
critical  concerns  when  multiple  arrays  operate  coherently.  A  five-subarray  DDSA  model 
is  used  to  examine  the  effects  of  fixed  errors  to  the  direction-of-arrivals  (DOAs) 
estimation  and  filling  method.  Each  subarray  is  comprised  of  10  elements  with  element 
spacings  equal  to  0A5A.  Subarray  center  distances  are  10 A.  Fixed  errors  are  uniformly 

o  o  o  o 

distributed  from  —21  to  21  (root  mean  square  (RMS)  values  from  0  to  12.1  ),  and  the 
SNR  is  varied  from  6  dB  to  21  dB  at  each  element.  Two  signals  with  equal  magnitude 

°  o 

and  non-coherent  phases  {  nl 5  and  -4;r/5)  from  DOAs  of  -10  and  15  relative  to 
broadside  are  impinging  on  the  DDSA. 

The  receiving  pattern  for  6  dB  SNR  for  the  ideal  contiguous  array  of  the  same 
aperture  size  as  DDSA,  original  DDSA  and  virtual  filling  method  are  compared  in  Figure 
5.  A  huge  improvement  in  terms  of  grating  lobes  and  sidelobe  suppression  can  be 
observed  after  applying  the  virtual  filling  method. 


Main  Beam  Scan  Angle  in  sin(o) 


Figure  5.  Pattern  comparison  of  contiguous,  original  DDSA  and  virtual  filling. 
Assuming  no  fixed  errors  and  with  6  dB  SNR  at  each  element.  Taylor  amplitude  taper 

( IT  =5,  SLL=— 30  dB)  has  been  applied. 
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In  order  to  quantify  the  effect  of  fixed  errors  on  the  DOA  estimations,  a  plot  that 
compares  the  root  mean  square  error  (RMSE)  of  the  DOA  versus  RMS  phase  error  for 
different  SNR  levels  is  shown  in  Figure  6. 


Figure  6.  RMSE  of  DOA  versus  RMS  phase  errors  from  0^  to  12.1°  for  different 

element  SNR. 


Many  of  the  sources  of  the  fixed  errors  can  be  compensated  for  by  pre-calculation 
or  pre-measurement,  but  there  will  still  be  some  residual  errors  after  correction.  We 

o 

consider  phase  errors  up  to  21  and  examine  how  they  degrade  the  radiation  pattern. 

o 

Figure  7  has  a  plot  of  the  pattern  of  the  worst  case  (21  fixed  error)  at  6  dB  SNR.  By 
comparing  Figure  7  with  Figure  5,  we  see  that  the  effect  of  the  fixed  error  on  the 
receiving  pattern  is  to  increase  the  side  lobe  level  and  lower  the  main  beam  by  0.6  dB  due 
to  the  increase  of  RMSE  of  the  DOA. 
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Figure  7.  Pattern  comparison  of  contiguous,  original  DDSA  and  virtual  filled  DDSA 

o 

for  21  fixed  error  and  with  6  dB  SNR  at  each  element.  Taylor  amplitude  taper  (JT  =  5, 

SLL=— 30  dB)  has  been  applied. 
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I.  INTRODUCTION 


The  phased  array  is  generally  the  antenna  architecture  of  choice  for  most  modern 
high-performance  radar  and  communication  systems.  Phased  arrays  consist  of  a 
collection  of  individual  antennas  that  are  geometrically  arranged  and  excited  (phased)  so 
as  to  provide  the  desired  radiation  characteristics. 

Traditionally  antenna  arrays  have  been  constructed  with  a  large  number  of 
radiating  elements  distributed  over  a  given  confined  surface  with  an  average  distance  to 
the  nearest  neighbor  no  larger  than  A/2 .  Advanced  Active  Electronically  Scanned  Array 
(AESA)  systems,  which  structurally  integrate  the  arrays  into  the  platform,  have  opened 
up  the  possibility  of  antenna  systems  consisting  of  separated  subarray  apertures  where  the 
aperture  is  split  into  two  or  more  subarrays  and  separated  by  a  relatively  large  distance. 
Within  each  subarray  the  average  distance  between  the  radiating  elements  is  still  not 
larger  than  A/2  but  the  distances  between  the  phase  centers  of  the  subarrays  are  much 
larger  than  that.  An  increase  in  angular  measurement  performance  can  potentially  be 
gained  by  increasing  the  measurement  base  by  dividing  the  aperture  into  two  or  more 
subarrays  and  pulling  them  apart.  This  should  improve  both  angular  accuracy  and 
resolution  [1,2]. 

Also,  when  multiple  systems  on  the  same  platform  are  using  arrays  with  increased 
bandwidths,  these  systems  can  share  arrays.  Another  factor  impacting  architecture  is  the 
platform  design  philosophy  for  military  applications,  which  has  changed  dramatically 
with  the  advent  of  stealth  technology  and  requires  reduced  platform  signatures.  Due  to 
the  stealth  requirement,  it  is  difficult  to  find  an  available  area  sufficient  for  a  large  array 
on  board  a  ship,  so  it  might  be  necessary  to  use  several  relatively  small  noncontiguous 
(separated)  areas  (subarrays)  and  then  process  the  received  signal  coherently. 

“Hastily  formed”  subarray  systems  are  reconfigurable  and  expandable  arrays  for 
emergencies  that  can  be  deployed  on  the  sides  of  hills,  buildings,  or  trucks.  Some  of  these 
applications  are  illustrated  in  Figure  1. 
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Subarrays  randomly  or 
aperiodically  disiributcd 
over  large  surface  areas 


Assemble  “ad  hoc"  arrays  for  air  traffic 
control,  air  defense,  SIGIN7 


Figure  1.  Assembled  “ad  hoc”  arrays  for  air  traffic  control,  air  defense,  signal 

intelligence  (SIGINT). 


A.  PERIODIC  DISTRIBUTED  SUBARRAYS 

Traditional  periodically  distributed  subarrays  (subarrays  whose  centers  are 
equally  spaced  )  form  a  long  baseline  and  are  capable  of  very  accurate  angular  location  of 
targets  [3,  4].  However,  collectively  combining  periodic  widely  separated  subarrays 
results  in  unacceptable  grating  lobes,  and  these  lobes  cannot  be  suppressed  using 
traditional  windowing  methods.  Grating  lobes  appear  in  the  visible  region  if  the  subarrays 
are  nearly  periodic,  and  they  are  unwanted  because  of  the  ambiguities  that  accompany 
them.  Besides  that,  the  interference  power  from  another  emitter  can  severely  degrade  the 
system  performance  when  the  emitter  is  in  the  direction  of  a  grating  lobe.  Similar 
degradation  can  occur  when  a  sector  of  clutter  return,  which  consists  of  reflections  of 
transmitted  power  from  obstacles,  is  received  through  the  grating  lobes.  The  depiction  of 
the  desired  signal,  interference,  clutter  and  receiver  noise  can  be  seen  in  Figure  2.  The 
signal-to-interference  ratio  (SIR),  the  signal-to-clutter  ratio  (SCR)  and  the  signal-to-noise 
ratio  (SNR)  are  defined  as  the  signal  power  over  the  interference  power,  the  signal  power 

over  the  clutter  power,  and  the  signal  power  over  the  noise  power,  respectively. 
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Interference  (emitter) 


Figure  2.  An  illustration  of  the  desired  signal,  interference,  clutter  and  receiver 

noise. 

It  is  possible  to  suppress  the  grating  lobes  through  control  of  the  subarray  factor. 
In  other  words,  the  nulls  of  the  subarray  pattern  are  placed  at  angles  where  the  grating 
lobes  occur.  However,  the  resulting  constraints  are  severe.  Some  additional  improvement 
is  achieved  if  separate  transmitting  and  receiving  antennas  are  used,  and  grating  lobes  are 
allowed  for  only  one  of  the  two  antennas.  The  main  limitation  is  that  the  configurations 
of  the  transmitting  and  receiving  antennas  need  to  be  different,  and  hence  more  space  is 
required  [5]. 

B.  RANDOM  AND  APERIODIC  DISTRIBUTED  SUBARRAYS 

Another  approach  to  reducing  grating  lobes  is  randomization  of  the  subarray 
positions.  A  contiguous  random  fractal  array  is  an  example  that  utilizes  randomness  to 
suppress  the  grating  lobes  [6].  An  aperiodic  array  is  somewhere  between  periodic  and 
random  arrays.  For  a  linear  array  with  unequal  element  spacing,  it  is  possible  to  mitigate 
large  grating  lobes  because  of  the  “space  tapering”  effect  from  the  non-uniform  element 
spacing  [7]. 
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c. 


DIGITAL  ANTENNAS 


Digital  antennas  are  widely  used  today  because  of  their  flexibility  in  signal 
processing  and  suitability  for  low  profile  array  design  compared  to  their  analog 
counterparts  [8].  An  example  of  digital  array  antenna  architecture  is  shown  in  Figure  3. 
The  important  feature  of  this  architecture  is  that  element  level  baseband  in-phase  (/)  and 
quadrature  ( Q )  data  are  available  to  a  central  array  processor  (digital  beamformer).  A 
digital  antenna  architecture  allows  for  a  wide  range  of  waveform  control  and  digital 
beamforming  on  both  the  transmitting  and  receiving  sides. 


Figure  3.  Generic  digital  antenna  architecture. 


When  subarrays  are  used  together  in  a  digital  antenna  system  there  are  several 
beamforming  options.  The  first  is  to  let  them  act  independently  (with  or  without  a 
common  time  and  phase  reference).  The  system  can  schedule  the  antenna  resources  to 
improve  the  performance  [9].  For  example,  a  radar  may  add  a  second  array  to  increase 
the  scan  volume  or  reduce  the  frame  time. 

A  second  approach  is  to  use  the  antennas  collectively  as  if  they  are  elements  of  a 
larger  array,  i.e.,  they  act  as  subarrays  in  a  large  array  [10].  This  arrangement  could  be 
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referred  to  this  as  a  “synthetic  array.”  Depending  on  the  performance  requirements  and 
subarray  locations  and  orientation,  not  all  subarrays  need  to  be  used  at  the  same  instance. 

D.  GEOMETRICAL  ARRANGEMENT 

The  geometrical  arrangement  and  individual  capabilities  of  the  subarrays  greatly 
affect  the  potential  performance  of  the  synthetic  array.  Fundamental  array  parameters 
include: 

•  The  number  of  subarrays. 

•  The  size  of  the  subarrays  and  their  shapes. 

•  Element  types  and  polarization. 

•  Scan  capability  (maximum  scan  angle,  scan  rates  and  planes  of  scan). 

•  Physical  arrangement  (the  distribution  of  elements  and  subarrays  in 

space). 

•  Linear,  planar,  volumetric. 

•  Periodic,  aperiodic  or  random. 

•  Rotation  and  tilt  relative  to  a  global  reference. 

Based  on  the  physical  arrangement,  periodic  array  structures  can  be  linear,  planar 
or  volumetric.  Linear  and  planar  structures  are  most  often  used  because  they  are 
relatively  easy  to  design,  manufacturing  and  integrate  into  platforms,  and  they  can 
provide  any  desirable  beamwidth  and  gain.  In  this  research  we  shall  assume  that  all  of  the 
arrays  are  planar,  but  they  can  be  tilted,  rotated  and  have  different  numbers  of  elements 
and  element  spacings.  An  example  of  distributed  linear  array  of  Ns  subarrays  in  the  x-y 
plane  is  shown  in  Figure  4.  We  consider  only  rectangular  subarrays;  the  extension  to 
other  shapes  (e.g.,  circular)  is  straightforward. 

Each  subarray  is  capable  of  operating  independently  (i.e.,  it  can  scan 
independently),  and  its  element  level  /  and  Q  data  are  available  to  a  central 
beamformer/processor.  A  primary  interest  in  this  research  is  to  develop  beamforming 
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techniques  on  the  receiving  side  because  there  is  more  flexibility  in  processing  the 
received  data.  The  flexibility  of  the  hardware  on  the  transmitting  side  is  generally  more 
limited. 


Airay  element  Subarray 


y 


m=  1  m=Ns 

x 

Figure  4.  Distributed  linear  array  structure. 

Periodic  spacing  of  array  elements  along  a  line  or  over  a  plane  has  been  a  major 
way  of  creating  a  phased  array  antenna  for  decades.  As  long  as  the  element  spacing  for 
the  array  is  less  than  a  half  wavelength,  grating  lobes  will  not  exist  in  the  far  field 
radiation  pattern  when  scanned.  An  advantage  of  using  periodic  structures  is  that  lobe  and 
null  locations  of  the  radiation  pattern  are  predictable,  and  they  result  in  lower 
manufacturing  and  replacement/maintenance  cost.  Also,  pattern  synthesizing  and 
beamforming  are  relative  easy  compared  to  that  of  non-periodic  structures. 

A  disadvantage  of  the  periodic  array  structure  is  that  grating  lobes  become  a 
problem  for  element  spacing  larger  than  a  half  wavelength.  Because  the  pattern  is 
determined  by  the  element  spacing  in  wavelengths,  narrow  bandwidth  is  not  uncommon. 
Generally  speaking,  the  aperture  size  must  be  increased,  and  more  radiation  elements 
(larger  area)  are  needed  to  decrease  the  beamwidth.  Larger  spacing  is  also  used  to  reduce 
the  number  of  elements  in  an  array  to  reduce  cost,  weight  and  volume. 

For  military  applications,  stealth  is  one  of  the  crucial  issues  that  cannot  be  ignored 
for  today’s  warfare.  A  revolutionary  idea  that  could  provide  the  potential  of  using  the 
entire  length  and  width  of  a  large  ship  as  the  aperture  is  to  let  the  stealth  designed  ship 
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structure  dictate  where  the  subarrays  can  be  placed.  The  aperture  can  be  used  collectively 
or  individually  for  various  purposes  such  as  communication,  surveillance,  guidance  and 
control,  and  even  threat  imaging  [5].  This  idea  of  having  integrated  subarray  elements  in 
the  ship’s  hull  forming  the  aperture  for  the  radar  is  not  entirely  new  [1 1],  but  there  arc 
problems  associated  with  a  distributed  array  system  (DAS)  such  as  those  shown  in  Figure 
1.  Not  only  are  there  grating  lobes  due  to  the  subarray  spacing,  but  the  distribution  of 
frequency  and  time  references,  and  calibration  and  error  correction  are  significant 
challenges. 

E.  DISSERTATION  OBJECTIVE 

Very  few  researchers  have  looked  at  using  a  large  number  of  arrays  collectively, 
in  particular,  if  the  subarrays  are  not  identical.  The  conventional  approaches  to  reducing 
grating  lobes  have  mainly  been  applied  to  individual  arrays  but  not  to  multiple  arrays 
with  gaps.  Likewise,  random  methods  are  often  applied  to  the  distribution  of  elements 
within  an  array  but  not  between  subarrays. 

A  primary  objective  of  this  research  is  to  investigate  the  wide  variety  of 
processing  techniques  that  can  be  applied  with  digital  distributed  subarray  antennas  for 
grating  lobe  suppression.  Many  radar,  communication  and  electronic  warfare  (EW) 
systems’  performance  and  capabilities  can  be  improved  by  using  multiple  antennas  (or 
equivalently  one  large  synthetic  antenna).  Potential  improvements  are  the  increase  in 
signal-to-noise  ratio,  high  angular  resolution,  wider  bandwidth,  and  reconfigurability. 

In  this  research,  we  focus  on  distributed  antennas  that  arc  comprised  of  subarrays 
that  can  operate  individually  or  collectively.  It  is  assumed  that  no  grating  lobes  appear  in 
the  visible  region  for  each  subarray  when  scanned.  We  develop  techniques  for  grating 
lobe  suppression  on  both  the  transmitting  and  receiving  sides  of  the  DAS.  Traditional 
solutions  and  new  methods  will  be  examined  in  detail  via  numerical  simulation  to 
quantify  the  performance  limitations  when  applied  in  combination.  One  contribution  of 
this  research  is  a  hybrid  approach  that  uses  a  combination  of  suppression  techniques  on 
both  the  transmitting  and  receiving  sides.  The  result  is  an  improved  two-way  pattern 
performance. 
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Another  contribution  is  the  development  of  new  receiving  processing  methods  to 
suppress  grating  lobes  and  improve  the  signal-to-clutter  ratio  and  signal-to-interference 
ratio.  A  final  contribution  shows  the  relationship  (and  hence  tradeoffs)  between  thermal 
noise,  array  errors,  and  the  grating  lobe  suppression  effectiveness.  The  consideration  of 
array  errors  addresses  the  issue  of  array  calibration  and  synchronization,  which  are 
critical  concerns  when  multiple  arrays  operate  coherently. 

F.  RELATED  WORK 

Steinberg  and  Yadin  in  1982  [12]  first  introduced  the  concept  of  a  distributed 
airborne  array  to  make  use  of  a  large  aperture  to  achieve  higher  angular  resolution.  Lin’s 
1983  paper  [3]  was  the  first  to  consider  using  coherently  internetted  mini-radars  to 
accomplish  surveillance  and  tracking  efficiently.  Lin’s  2003  thesis  [5]  adapted  the 
concept  of  distributed  subarray  antennas  and  proposed  it  for  both  Multi-function  Array 
Radar  (MFAR)  and  Very  High  Frequency  (VHF)  applications.  By  combining  distributed 
subarrays  on  the  available  areas  of  a  constrained  platform,  the  MFAR  or  VHF  distributed 
subarray  antennas  (DSAs)  can  achieve  the  maximum  resolution  and  potential  reductions 
in  cost  and  complexity. 

Many  methods  have  been  employed  to  reduce  or  eliminate  the  grating  lobes  for 
DSAs,  but  all  have  their  limitations  and  disadvantages.  (We  refer  to  these  collectively  as 
“traditional”  or  “conventional”  techniques.)  Within  a  single  array,  a  common  approach  is 
to  place  subarray  nulls  at  grating  lobe  locations  using  overlapping  subarrays  [13],  but  this 
severely  limits  the  array  geometry.  Another  approach  is  to  rotate  or  tilt  the  subarrays, 
thereby  reducing  the  periodicity  [14].  The  grating  lobe  level  varies  as  201og(l/V  ) ,  where 

Ns  is  the  number  of  subarrays.  To  be  effective,  this  method  requires  a  large  number  of 
subarrays.  Random  or  fractal  element  spacings  within  the  subarrays  and  randomizing  the 
number  of  elements  between  subarrays  have  been  used  [15].  Again,  large  numbers  of 
elements  and  subarrays  are  needed  for  truly  random  behavior,  and  only  modest  grating 
lobe  suppression  is  achieved  [16].  (For  a  128  element  linear  array,  the  improvement  is 
about  6  dB.)  Multiplicative  beamforming  has  also  been  applied  to  suppress  grating 
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lobes  [17],  but  the  resultant  loss  in  SNR-gain  (an  average  of  6  dB)  is  the  main  drawback 
of  this  method  [18]. 

All  of  these  methods  adopts  only  a  single  solution  for  grating  lobe  suppression, 
and  their  effectiveness  is  generally  not  good.  An  alternative  approach  to  improving  the 
grating  lobe  suppression  is  a  hybrid  method,  which  includes  both  a  physical  treatment  in 
tandem  with  digital  signal  processing  (DSP).  On  the  transmitting  side  ‘'mild” 
randomization  of  subarrays  is  used  to  lower  the  grating  lobes.  On  the  receiving  side,  we 
propose  virtual  filling  of  the  gaps  between  the  distributed  subarrays  to  eliminate  the 
grating  lobes  on  the  receiving  side  so  that  the  response  of  a  single  large  contiguous  array 
is  synthesized.  This  approach  was  recently  suggested  in  [19]  to  fill  gaps  in  the  array 
matrices  for  super-resolution  direction-of-arrival  (DOAs)  estimation.  They  use  minimum 
weighted  norm  (MWN)  and  super  spatially  variant  apodization  (Super-SVA)  for  virtual 
filling,  which  requires  significant  computational  power  and  large  numbers  of  time 
snapshots  (time  samples).  Super-SVA  performance  degradation  occurs  for  coherent 
signals. 

For  the  virtual  filling  method,  DOA  estimation  is  crucial  to  synthesizing  the 
virtual  element  weights.  Super  resolution  techniques  that  are  based  on  the  eigen-structure 
of  the  input  covariance  matrix,  such  as  multiple  signal  classification  (MUSIC)  [20],  root- 
MUSIC  [21]  and  estimation  of  signal  parameters  via  rotational  invariance  techniques 
(ESPRIT)  [22],  are  used  to  generate  the  high  resolution  DOA  estimates.  However,  an 
extra  step  of  spatial  smoothing  is  needed  for  correlated  or  coherent  signals,  and  a  large 
number  of  snapshots  are  required  for  acceptable  accuracy.  Sarkar  [23]  utilized  the  matrix 
pencil  (MP)  method  to  get  the  DOA  of  the  signals  in  a  coherent  multi-path  environment. 
In  the  MP  method,  based  on  the  spatial  samples  of  the  data,  the  analysis  is  done  on  each 
snapshot,  and  therefore  non-stationary  environments  can  be  handled  easily.  In  recent 
years,  the  MP  method  has  received  more  attention  due  to  its  computational  simplicity  and 
accuracy  in  noisy  environments.  Comparative  studies  with  other  DOA  estimation 
algorithms  can  be  found  [24-32]. 

However,  traditional  MP  is  designed  for  the  single  periodic  array  and  cannot  be 

used  directly  for  DSA  configurations.  Therefore,  to  meet  the  needs  of  this  research,  we 
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have  developed  a  modified  MP  method  for  DSA  and  devised  an  enhanced  three-step 
modified  MP  method  for  further  improvement  in  the  resolution  of  DOA  estimation. 

G.  ORGANIZATION  OF  THE  DISSERTATION 

First  a  survey  of  fundamental  grating  lobe  reduction  approaches  commonly  used 
on  single  sparse  arrays  (i.e.,  the  “conventional”  or  “traditional”  approaches)  is  presented. 
Then  these  methods  arc  extended  and  applied  to  the  DSA  problem  with  numerical 
simulations  to  assess  their  effectiveness.  Next  the  virtual  filling  solution  is  applied  on  the 
receiving  side,  and  finally  the  impact  on  the  two-way  (transmitting  and  receiving)  pattern 
is  examined. 

An  introduction  to  general  formulation  for  the  pattern  of  a  DSA  is  discussed  in 
Chapter  II.  Because  the  radiation  pattern  is  primarily  impacted  by  the  physical 
arrangement  of  radiation  elements,  i.e.,  array  lattices,  some  special  cases  are  discussed. 
They  include  planar  arrays  of  identical  subarrays  in  a  rectangular  grid  with  and  without 
random  displacements. 

There  are  a  number  of  traditional  approaches  to  lowering  the  grating  lobes  in  a 
single  sparse  periodic  array,  and  these  are  discussed  in  Chapter  III.  Also  included  in  this 
chapter  is  an  overview  of  combined  solution  approaches  and  the  various  advantages  and 
disadvantages  of  these  traditional  solutions. 

Grating  lobe  suppression  with  traditional  methods  is  considered  in  Chapter  IV. 
The  physical  arrangement  of  subarrays  and  the  multiplicative  beamforming  method  arc 
discussed  in  detail,  and  simulation  results  are  provided  for  each.  Because  of  the  limited 
performance  of  each  method,  a  hybrid  approach  which  combines  several  different 
traditional  methods  is  proposed  and  evaluated  in  this  chapter. 

The  new  filling  method  for  receiving  beamforming  is  considered  in  Chapter  V. 
The  DOA  estimation  technique  for  lion-stationary  environments  is  studied  and  tailored 
for  the  DSA  models.  This  is  done  first  from  the  classical  MP  method.  Next  MP  is 
modified  to  be  used  for  DSA  problems.  These  MP  studies  comprise  the  basis  of  the 
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proposed  filling  methods  for  grating  lobe  suppression,  even  though  they  could  stand 
alone  as  a  separate  research  problem. 

Potential  applications  of  proposed  methods  arc  discussed  in  Chapter  VI.  An 
extended  three-step  MP  method,  which  improves  the  estimation  accuracy,  is  presented. 
To  complete  a  system-wide  study,  both  the  transmitting  and  receiving  sides  must  be 
considered.  By  combining  random  subarray  sizes  on  the  transmitting  side  and  the  virtual 
filling  method  on  the  receiving  side,  a  novel  two-way  DSA  pattern  with  suppressed 
grating  lobes  is  achieved. 
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II.  THEORETICAL  BACKGROUND 


A.  A  GENERAL  FORMULATION  FOR  THE  PATTERN  OF  DISTRIBUTED 
DIGITAL  SUBARRAYS 

To  begin  we  present  the  general  formulation  of  radiation  patterns  for  distributed 
digital  subarrays  (DDSAs).  These  formulas  apply  to  any  arbitrary  three-dimensional 
arrangement  of  arrays.  Formulas  for  special  cases  such  as  planar  and  random  subarrays 
will  also  be  covered  in  this  chapter.  The  spherical  coordinate  system  that  is  used 
throughout  is  shown  in  Figure  5. 


z 


Figure  5.  Spherical  coordinate  system  for  DDSA. 
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1. 


General  Formulas 


In  Figure  6,  an  illustration  of  the  general  array  geometry  is  shown.  The  subarrays 
are  defined  by: 

Ns  =  number  of  subarrays, 
m— subarray  index,  m  =  1, 2, and 

xs{m)9  ys{m)9  zs(m)  =  coordinates  of  subarray  m  in  the  global  system. 


The  subarrays  can  be  rotated  and  tilted  with  respect  to  the  global  origin.  In  the  global 
system  (#,  (/))  the  direction  cosines  are 

w  =  sin#cos^ 

v  =  sin#sin^  (1) 

w  =  cos  6. 

The  scan  direction  in  global  coordinates  is  ( 0S9  /s)9  with  direction  cosines 

us  =  sin  0S  cos  (fs 

vs  =  sin  0S  sin  <j)s  (2) 

ws  -  cos0s. 

In  the  local  subarray  m  coordinate  system  (0m,  /m)  the  direction  cosines  are 

u =  sin  em  cos 

vm  -  sin  0m  sin 4  (3) 

wm  ~  COS  0m , 

and  the  scan  angle  direction  cosines  in  the  local  subarray  coordinate  system  are 

usm  =  sin  0sm  C°s4„ 

vsm  =  sin  9sm  sin  </>sm  (4) 

Wsm=  COS  0sm. 

A  rotation  matrix  can  be  used  to  obtain  the  subarray  direction  cosines  from  the  global 
ones,  and  vice  versa. 
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Figure  6.  Subarray  m  and  its  local  coordinate  system  relative  to  the  global  origin. 

Consider  planar  rectangular  subarrays  with  all  arrays  in  the  local  z  =  0  plane.  The 
element  spacing  of  the  subarrays  can  be  different,  but  the  spacing  within  each  subarray  is 
equal.  The  number  of  elements  in  the  subarrays  can  vary.  Define: 

Nx  (m),  N y  (m )  =  number  of  elements  in  the  local  x  and  y  directions  for  subarray  m, 

N(m)  =  Nx{m)N y(m)  =  total  number  of  elements  in  subarray  w, 

dx(m),dy(m)  =  spacing  between  elements  for  subarray  m, 

Lx{m )  =  Nx(m)dx(m)  =  length  of  the  mth  subarray  in  its  x  direction,  and 

Ly(rn)  =  Ny(m)dy{m)  =  length  of  the  mlh  subarray  in  its  y  direction. 

We  also  define  k  =  2k IX ,  where  X  is  the  wavelength  at  the  frequency  of  operation. 

The  complex  pattern  of  the  mlh  subarray  is  given  by  the  sum  (  e,cor  time 
dependence  assumed  and  suppressed) 

Nx(m)Ny(m)  . 

Fs(m)=  Z  I  a(rn,p,q)ew{m'p'qUxp[jk(x(m,p)um  +y(m,q)vm)]  (5) 

P= 1  ^=1 
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where 


w(m9p9q)  =  a(m9  p9q)e^^m,p^  (6) 

is  the  complex  weight  (amplitude  a  and  phase  y/  )  at  element  p ,  q  of  subarray  m  applied 
for  scanning,  side  lobe  control,  beam  shaping  and  error  compensation.  For  equally  spaced 
elements,  with  each  subarray  centered  at  its  local  origin 

2  P  ( A^*W)  +  1)  dx  (m)  =  P(m,  p)dx  (m) 

(7) 


x(m,  p)  =  ■ 


2q-(N  v(m)  +  \) 

y(m,  q)  = - - - dy  (m)  =  Q(m,  q)d y  (m) 


so  that 


Nx(m)Ny(m) 

Fs(m)=  I  Z  w(m,p,q)exp[jk(P(m,p)dx(m)um  +Q(nt,q)dy(m)vm)j.  (8) 
p= i  q= i 

Note  that  when  receiving,  the  exponential  factor  would  be  obtained  from  the  element 
baseband  /  and  Q  samples: 

I(m,p,q)  +  jQ(m,p,q)  =  exp[jk(P(m,p)dx(m)um  +Q(m,q)dy{m)vm)\.  (9) 
The  array  factor  (constructed  from  the  centers  of  each  subarray)  is 


Fa=  Z  exp\_jk(xs(m)u  +  ys(m)v  +  zs(m)w)~\. 
m= 1 

The  total  distributed  array  factor  is  given  by 

Ns 

F  =  £  F9(/i7)exp[yA:(x9(m)i/  +  79(m)v  +  z9(m)w)] 

m= 1 

Ns  [ Nx(m)Ny(m ) 


(10) 


=  Z 


m— 1  I  p=\  q= 1 


Z  \v(m9p,q)exp[jk(P(m,p)dx(m)um  +  Q(m,q)d l  (11) 


x  exp£  jk  ( xs(m)u  +  ys{m)v  +  zs(w)w)J. 


If  the  weights  are  separable  in  the  x  and  y  coordinates  such  that 

w(m,p,q)  =  wx(m,p)wy(m,q), 


(12) 


then 
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Ns  [*,(*)  .  ,Ny(m)  r 

F=  I  j  S  wx(m,p)cxp[jkP(m,p)dx(m)um\  I  Wj, (m, <?) exp[y/cQ(m, q)dy(rn)vm J 
m= 1  [  p= 1  g=l 

x  exp  [  jk  ( (m)w + ys  (m)v  +  z5  (m)w)] . 

(13) 

Furthermore,  if  the  amplitude  weights  are  uniform  (wx  =  wy=  l)and  a  linear  phase  is 
added  to  scan  the  beam  in  the  local  coordinate  system  (usm9  vsm) ,  then 


x  expj^  jk  [xs{m)u  +  ys(m)v  +  zs(m)w)_ |. 


In  order  to  linearize  the  phase  across  the  entire  array  a  scanning  phase  (wv,  yv,  ws) 
should  be  added  to  the  exponential  factor  in  Eq.  (14): 

exP  {jk  [ xs(m)(u  -  us)  +  ys(m)(v  -  vs)  +  zs(in)(iv-ws)]}.  (15) 

However,  these  scanning  phases  could  be  added  to  the  element  weights  as  constant  values 
for  each  subarray. 

To  complete  the  expression  for  the  pattern  an  element  factor  must  be  added.  In  the 
local  subarray  coordinates  ( 0m ,  (f)m)  the  element  factor  S  for  subarray  m  can  be  expressed 
as 

S(m,em,<j>m)  =  Se(m,0m ,  <fm  )0m  +  (m,  0m,fm  )j>m .  (16) 

This  allows  for  the  possibility  that  elements  are  different  for  each  subarray,  which 
generally  would  not  be  the  case.  The  final,  most  general  expression  for  the  total  array 
pattern  is 

ffs 

F(u,v)=  I  Fs(m)  S(rn,um,vm) 

m= 1  (17) 

x  exp  {jk  [x9(w)(w  -  us)  +  ys(m)(v  -  v?)  +  zs(m){w-  wv)]}. 
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The  normalized  power  pattern  is  computed  by 


Jnonn(«»v)  = 


F(u,v)  F(u,v)»  F(u,v)* 


F{u,v )2 

max 


F(u,v) 


(18) 


max 


The  array’s  gain  pattern  can  be  expressed  as  the  peak  gain  times  this  normalized  power 
pattern 

G(u,v)  =  G0Pnorm(u,v)  (19) 

where  Gq  is  the  main  beam  gain. 

2.  Special  Case  1:  Planar  Array  of  Identical  Subarrays  in  a  Rectangular 
Grid 

a.  Specialized  Formulas 

The  index  m  can  be  dropped  from  the  subarray  quantities  if  all  subarrays 
are  identical.  As  shown  in  Figure  7,  let:  Nsx ,  Nsv  =  number  of  subarrays  in  the  x  and  y- 

directions  (Ns  =  NsxNSy),  and  sXo,Sy0  =  center-to-center  spacing  between  subarrays  in 

the  x  and  y  directions.  Then 

-2P-{Nsx+iy 


xs(p)  = 
ys(q) 


2 

2<y-(A',,+i) 


•%  =  p(P)sx0 

syo=Q(q)syo. 


(20) 


The  subscript  m  can  also  be  dropped  from  the  direction  cosines  because  all 
subarrays  are  aligned  with  the  global  coordinate  system.  Now  the  array  factor  can  be 
reduced  to  closed  form 


Nsx  NSy 


Fa  =  f  Z  expry*(P(p)jjKb(ii-i/J)  +  e(^)5J,0(v-v(S))l 

p=\q=\ 


sin 


N 


sx 


ksXo(u-Us) 


sin 


N, 


sy 


ksyAv~vs ) 


2  y° 


(21) 


ksY  , 

ks„ 

sin 

2°(m  us) 

sin 

o 

< 

1 

r*^ 
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Figure  7.  Periodic  array  in  a  rectangular  grid. 


The  total  complex  pattern  is  constructed  from  the  product  of  the  subarray 
factor  and  array  factor: 


b.  Lattice  Grating  Lobes 

The  situation  of  interest  is  one  where  the  element  spacing  is  small  enough 
so  that  individual  subarray  grating  lobes  do  not  occur  for  any  scan  angle.  Any  grating 
lobes  that  occur  are  due  to  the  gaps  between  the  subarrays  in  the  rectangular  grid  (i.e.,  the 
lattice  grating  lobes,  also  called  construction  grating  lobes).  Several  conventional 
techniques  have  been  used  to  suppress  the  lattice  grating  lobes: 

•  Employ  an  array  element  that  has  nulls  at  the  locations  of  the 
grating  lobes:  This  approach  is  limited  to  a  couple  of  wide  angle 
grating  lobes  at  best. 

•  Employ  other  periodic  lattice  configurations:  For  example,  a 
triangular  lattice  shifts  the  grating  lobes  off  of  the  principal  planes 

19 


(in  direction  cosine  space).  Circular,  hexagonal,  octagonal,  etc. 
grids  have  been  used.  Other  shapes  have  also  been  investigated, 
such  as  fractal  and  spiral  constructions. 

•  Adjust  the  subarray  size  and  element  spacing:  Select  the  subarray 
spacing  and  element  spacing  so  that  subarray  pattern  nulls  fall  at 
the  grating  lobe  locations.  Again,  there  are  limitations  due  to  the 
constraints  on  the  size  of  the  subarrays  and  the  gaps  between  them. 

•  Multiplicative  beamforming:  A  single  “master”  subarray  pattern  is 
used  to  suppress  grating  lobes,  and  hence  it  reduces  the  design 
constraints.  However,  it  has  other  implications  on  the  beamforming 
performance  that  will  be  discussed  in  Chapter  II,  Section  B.  It  can 
only  be  applied  on  the  receiving  side. 

•  Minor  perturbation  of  the  geometry:  Minor  changes  in  the 
geometry  such  as  rotation  and  tilt  of  the  subarrays  can  be  used  to 
reduce  (but  not  eliminate)  the  grating  lobes.  The  perturbations 
complicate  the  design,  deployment  and  manufacturing  of  the  array. 
The  changes  also  reduce  gain  and  increase  average  side  lobe  levels, 
but  these  consequences  are  usually  tolerable. 

•  Aperiodic  or  random  array  configurations:  If  the  array  has  no 
periodicity  there  will  be  no  lattice  grating  lobes.  Taken  to  the 
extreme  this  would  be  a  random  array,  where  the  subarray 
locations  are  completely  random.  Most  applications  do  not  allow 
for  randomly  distributed  locations.  One  exception  would  be  a 
“randomly  thinned”  array  because  it  is  deployed  over  a  well- 
defined  area. 
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3.  Special  Case  2:  Planar  Array  of  Identical  Subarrays  in  a  Rectangular 
Grid  with  Random  Displacements 

In  this  case  xs{ni),ys{m)  are  random  numbers  with  a  known  probability  density 
function  (PDF),  and  hence  a  known  average  and  variance.  We  will  assume  there  is  no 
displacement  in  the  z  direction,  but  a  random  deviation  in  x  and  y  can  occur  for  each 
subarray.  It  is  assumed  that  the  errors  are  small  enough  so  that  there  is  no  overlap  in  the 
subarrays  in  their  displaced  positions.  The  spacing  between  subarrays  can  be  written  as  a 
sum  of  the  error  free  distance  plus  the  random  error  in  the  x  and  y  directions  8x ,  Sv , 

respectively.  The  PDF  of  Sx ,  Sv  can  be  inferred  by  the  method  in  which  the  subarrays 

are  distributed.  These  are  generally  zero  mean  with  variances  S 2  and  S2  . 


The  location  of  the  subarray  p ,  q  is: 

xs(p,q)  =  2p  yv-v  +  -VA-0  +  Sx(p,q)  =  P(p)sXq  +  Sx(p,q) 

1  (23) 

2q-{N,v  +  \) 

ys(P’<l)  = - — - Sy0  +  Sy{p,q)  =  Q(q)syo  +  Sy{p,q). 

Note  that  xs ,  ys  are  now  functions  of  both  indices  p  and  q.  This  is  because  the  random  x 
and  y  displacements  must  be  assigned  independently  for  each  subarray  (i.e.,  not  just  for 
rows  and  columns)  so  that  the  errors  at  subarrays  are  uncorrclated. 


Now  we  compute  the  mean  power  pattern.  To  simplify  the  equations  we  consider 
only  the  x  direction  (which  would  be  the  case  for  a  linear  array  along  x).  The  result  can 
easily  be  extended  to  the  y  direction  by  analogy.  The  mean  power  pattern  is  the  expected 
value  of  the  power  pattern: 


F(u)*F(u)*=  f  f  exp T jk ( P(p)sXo  +  Sp )(u-us)  exp[ -jk(p(n)sXo+S„)(u-us 

P= 1  77  =  1  L  J  L 


N'. 


=  f  f  exp[y'A:5ro  (P(  p)  -  P{n))(u  -  us )]  exp  [/A-  ( Sp  -  Sn )  (;/  -  us )]. 

P=\  77  =  1 


(24) 
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Wc  can  define  a  new  random  variable  A  =  k{Sp  -  Sn)(u  -  us)  .  The  mean  and 

variance  of  this  new  random  variable  (0,  A2 )  can  be  determined  from  the  mean  and 
variance  of  S  and  5 

p  n 

A2  =  k2  [S2p  +  S2^(a-us )' .  (25) 

When  p=n  in  Eq.  (24)  all  exponents  are  zero  and 

Z  Z  (1)  =  Nsx.  (26) 

P  =  ” 

Also,  using  the  fact  that  the  maximum  value  of  F  will  be  Nsc 

lF(w)|2norm"T72“|iV-«+  ^  t  cxp\jk  sXQ(P(p)  -  P(n))(u  -  us)]elA  1.  (27) 

Nsx  !  P* «  J 


The  expected  value  of  the  exponential  is  [33] 


p  =  n 
p  ^  n. 


(28) 


Now,  by  adding  and  subtracting 


N„ 


t  Z(l) 


p= " 


(29) 


and  using  the  +  with  the  second  term  in  the  curly  brackets  in  Eq.  (27)  to  complete  the 
sum,  we  get 


norm 


(30) 


where  Fa^  is  the  normalized  error  free  (unperturbed)  pattern.  As  the  error  increases,  the 

energy  in  the  error  free  pattern  is  transferred  to  the  second  term,  which  represents  random 
side  lobe  ‘"noise.” 


22 


Equation  (30)  extended  to  include  they  dimension  for  a  planar  array  is 
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Fa(u,vY 


norm 


Ns  -  e~^2  Ns  +  e~^ 


FaSu,v) 


=  e 


-A 


Fa0(u’v) 


+ 


1-e 


-A* 


(31) 


Note  that  the  subarray  and  element  factors  still  need  to  be  included  to  get  the  total  array 
pattern. 


B.  PERIODIC  DISTRIBUTED  DIGITAL  SUBARRAYS 

1.  Background 

As  mentioned  in  Chapter  I,  more  arrays  are  used  on  the  same  platform  to  facilitate 
the  multi-function  requirements  of  today’s  warfare.  Those  arrays  might  be  distributed 
over  the  platform  in  the  form  of  periodic  spacing.  Each  array  itself  is  an  element  of  a 
periodic  array  (i.e.,  a  subarray),  and  depending  on  the  mission  type,  subarrays  can  be 
operated  independently  or  collectively.  However,  collectively  operated  distributed 
subarrays  provide  a  narrower  main  beam  beamwidth,  hence  higher  angular  resolution  and 
better  SNR,  compared  to  the  sum  of  individually  operated  subarrays. 

The  pattern  function  for  the  case  of  a  rectangular  grid  of  subarrays  was  presented 
in  Eq.  (22).  It  was  noted  that  a  large  spacing  between  subarray  centers  leads  to  a  major 
disadvantage  of  this  kind  of  arrangement,  which  is  grating  lobes. 

In  the  following  section,  two  methods  are  proposed  to  eliminate  or  suppress  the 
grating  lobes  that  appear  in  the  periodic  distributed  subarrays.  The  first  method, 
multiplicative  beamforming,  has  been  used  in  radio  astronomy,  dealing  with  subarrays 
with  large  element  spacing  in  one  of  the  two  subarrays.  The  second  method,  subarray 
rotation,  has  been  applied  to  contiguous  arrays  for  grating  lobe  suppression. 

2.  Multiplicative  Beamforming  Applied  to  Distributed  Subarrays 

The  multiplicative  beamforming  method  has  been  primarily  used  on  thinned 
arrays  in  order  to  use  fewer  antenna  elements  yet  achieve  high  resolution.  A  major 
advance  in  this  method  was  made  by  Davies  and  Ward  [34],  who  introduced 
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interdependent  amplitude  tapers  for  two  constituent  subarrays  to  synthesize  a  desired  low 
side  lobe  radiation  pattern.  The  thinning  rates  that  can  be  achieved  are  of  the  order  of 
80%.  The  performance  of  multiplicative  receiving  systems  in  real  digital  beamforming 
(DBF)  radar  measurements  was  studied  in  [17].  A  compound  interferometer  is  a 
multiplicative  receiving  array.  Its  product  pattern  can  be  made  equal  to  the  power  pattern 
of  a  conventional  linear  array  with  uniform  element  weighting  was  shown  in  [35]. 
MacPhie  [36]  demonstrated  the  Mills  Cross  can  be  used  to  obtain  a  multiplicative  pattern 
equal  to  the  power  pattern  of  a  conventional  planar  array  but  with  fewer  elements.  An 
improvement  of  the  azimuthal  resolution  of  high-frequency  surface-wave  radar  (HFSWR) 
with  this  method  has  been  demonstrated  in  [37].  Lili  Xu  et  al.  [38]  utilized  the  concept  of 
the  multiplicative  array  to  improve  the  performance  of  nulling  antennas  as  well. 

The  primary  idea  of  multiplicative  beamforming  is  to  use  two  subarrays,  one  with 
element  spacing  of  /1/2  and  the  other  a  thinned  array  with  element  spacing  larger  than 
A/2  as  shown  in  Figure  8.  The  composite  directional  response  of  this  thinned 
multiplicative  array  can  be  seen  in  Figure  9.  These  subarrays  are  built  such  that  the  zeros 
of  the  first  subarray  match  the  grating  lobes  of  the  second  subarray.  This  processing 
implies  the  loss  of  SNR  gain,  defined  as  the  ratio  of  the  output  SNR  over  the  input  SNR, 
as  the  energy  in  the  grating  lobes  is  canceled.  It  is  only  used  on  the  receiving  side  because 
of  the  multiplication  requirement  [18,  37]. 


Figure  8.  Thinned  multiplicative  array  configuration  with  coincident  subarrays 

(from  [34]). 
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directional  response  caused 


Figure  9.  Composite  directional  responses  from  a  thinned  multiplicative  array  with 

coincident  subarrays  (from  [34]). 

The  other  limitation  of  this  method  is  the  problem  of  multiple  target  cross 
products  which  tend  to  provide  confusing  directional  information  for  complex  multi- 
target  excitation  of  the  array.  However,  if  the  signal  sources  are  uncorrelated,  this  cross 
product  averages  to  zero.  For  radar  applications,  the  adoption  of  frequency  agile 
transmissions  or  Doppler  processing  can  serve  to  decorrelate  the  returns  so  that  the  cross 
products  can  be  reduced  by  integration  [34]. 

3.  Distributed  Subarray  with  Subarray  Rotation 

Non-uniform  spacing  between  elements  would  remove  grating  lobes  but  normally 
at  the  cost  of  increased  side  lobes.  Optimization  methods  that  try  to  reduce  the  side  lobes 
require  substantial  computation  time  and  power  [39].  Random  designs  become  practical 
only  when  the  number  of  elements  is  very  large,  typically  1000  or  more.  A  new  method 
proposed  by  Agrawal  [14]  can  be  used  to  suppress  the  grating  lobes  to  a  certain  extent  by 
simply  rotating  subarrays  with  respect  to  each  other  by  specified  angle.  The  subarray 
arrangement  and  the  simulation  results  are  shown  in  Figure  10  and  Figure  11, 
respectively. 
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o  o  o  o 

Figure  10.  Eighty-element  array  with  four  subarrays  rotated  by  0  ,15  ,30  and  15 

(from  [14]). 


Figure  1 1.  Computed  radiation  patterns:  (a)  Eighty-element  array  without  subarray 
rotation,  (b)  Eighty-element  array  of  Figure  10  (  cl  =  2.6/1)  (from  [14]). 
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In  this  example,  the  array  with  a  triangular  lattice  is  divided  into  four  subarrays. 
The  element  spacing  d  is  equal  to  2.6/1  and  hence  grating  lobes  exist.  By  physically 

o  o  o  o 

rotating  subarrays  by  0  ,  15  ,  30  and  -15  ,  grating  lobes,  which  remain  at  the  same 
angular  distance  from  the  main  beam,  multiply  in  number  by  the  number  of  subarrays 
while  their  amplitude  is  divided  by  the  same  number  [14].  Kerby  [40]  examined  the 
behaviors  of  the  periodic  array  of  random  subarrays,  arrays  of  periodically  rotated 
random  subarrays,  and  arrays  of  randomly  rotated  random  subarrays  and  demonstrated 
the  lowered  side  lobe  level  of  the  array  factor  by  subarray  rotation. 

This  approach  can  be  extended  to  periodic  distributed  subarrays  to  lower  the 
magnitude  of  the  grating  lobes. 

C.  APERIODIC  AND  RANDOM  DISTRIBUTED  DIGITAL  SUBARRAYS 

1.  Background 

Aperiodic  and  random  arrays  provide  a  possibility  of  using  fewer  elements  than  a 
periodic  array  to  achieve  the  required  side  lobe  level  (SLL),  bcamwidth  (BW)  and 
directive  gain.  By  breaking  the  periodicity,  grating  lobes  can  be  mitigated  or  suppressed. 
Randomly  thinning  an  array  is  one  method  used  to  accomplish  this.  For  example,  every 
element  in  a  filled  array  can  be  assigned  a  number  between  0  and  1  from  a  uniform  PDF. 
If  the  array  is  to  be  thinned  by  50%  then  elements  with  numbers  less  than  0.5  could  be 
removed. 

Antenna  arrays  with  randomly  spaced  elements  have  been  studied  in  [41].  The 
probabilistic  properties  of  an  antenna  array  when  its  elements  are  placed  at  random  over 
an  aperture  according  to  a  given  distribution  were  presented.  It  was  found  that  the 
required  number  of  elements  is  closely  related  to  the  desired  side  lobe  level  and  is  almost 
independent  of  the  aperture  dimension.  The  resolution  or  the  bcamwidth  depends  mainly 
on  the  aperture  size.  The  directive  gain  is  proportional  to  the  number  of  elements  used  if 
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the  average  spacing  is  large.  As  a  consequence  the  number  of  elements  required  is 
considerably  less  than  one  with  uniform  spacing. 

In  practice,  taking  the  manufacturing  deviations  and  element  position  errors  due 
to  external  forces  into  account,  no  subarrays  are  identical.  Those  “errors”  are  usually 
random  variables  and  can  only  be  characterized  in  terms  of  probabilities  or  distributions. 
Thus,  in  a  sense,  all  practical  arrays  have  some  randomness,  but  the  errors  are  controlled 
so  that  they  are  small  compared  to  the  wavelength. 

Another  advantage  of  random  phased  arrays  is  the  potential  for  bandwidth 
improvement,  i.e.,  wideband  operation.  Aperiodicity  is  generally  recognized  as  one 
effective  way  to  extend  the  useful  bandwidth  of  antenna  arrays.  Random  arrays  can  be 
designed  with  very  little  pattern  variation  and  no  grating  lobes  over  wide  frequency 
ranges  [42].  Goffer,  et  al.  [43]  discussed  the  implementation  of  divided  arrays  with 
random  sizes  and  random  locations  of  their  centers  in  order  to  avoid  periodicity.  Closed- 
form  expressions  for  the  average  array  factor,  the  average  power  pattern  and  its  variance 
are  presented. 

However,  lack  of  periodicity  makes  manufacturing  more  complex  and  costly. 
Maintenance  will  be  another  issue  if  truly  random  arrays  are  used.  A  periodic  array  with 
random  subarrays  has  been  proposed  as  a  way  to  reduce  the  manufacturing  cost  but  retain 
some  properties  of  random  array  [40].  We  would  like  to  further  lower  the  manufacturing 
cost  and  increase  the  functionality  of  each  subarray  by  using  periodic  subarrays  with 
random  location  perturbations  or  periodic  subarrays  with  random  sizes. 

2.  Random  Array 

Collaborative  beamforming  for  wireless  sensor  networks  (WSNs)  has  been  used 
to  increase  the  transmission  range  of  individual  sensor  nodes  [44].  A  random  array  of 
subarrays  can  be  considered  as  a  wired  sensor  network.  The  individual  signals  from 
sensor  nodes  arrive  coherently  and  add  constructively  in  the  intended  direction.  The 
radiated  power  is  concentrated  in  a  certain  direction  and  reduces  the  power  in  other 
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directions.  The  beam  pattern  characteristics  of  collaborative  beamforming  have  been 
recently  derived  in  [45]  using  the  random  array  theory  and  assuming  that  sensor  nodes  in 
one  cluster  of  the  WSN  are  uniformly  distributed  over  a  circular  disk. 


Another  paper  derived  the  beam  pattern  for  the  WSN  with  Gaussian  distributed 
sensor  nodes  [46].  The  average  power  pattern  for  N  uniformly  distributed  nodes 
is  [45,46]: 
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where  Jn(x)  is  the  nlh  Bessel  function  of  the  first  kind.  The  radius  of  the  disk 

normalized  by  the  wavelength  is  a((ff)  =  4/rftsin(^/2)  and  R  =  R/A.  For  Gaussian 
distributed  nodes: 
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where  cc  =  4;rsin(^/2)  and  (j2  =  (jII/ l2  is  the  normalized  variance  of  the  Gaussian 
distribution. 


The  average  power  pattern  of  a  uniformly  distributed  WSN  is  shown  in  Figure  12. 
Several  values  of  R  with  N  =  16  and  256  are  plotted  for  comparison.  As  can  be 
observed,  the  side  lobe  approaches  1  IN  as  the  beam  angle  moves  away  from  the  main 
beam  as  predicted  by  Eq.  (32).  As  N  increases,  the  peak  side  lobe  level  goes  down.  The 
main  beam  width  is  inversely  proportional  to  R  . 
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Figure  12.  Average  beam  pattern  with  different  R  and  N  =  16  and 

256  (from  [45]). 


3.  Random  Distributed  Subarrays 

a.  Periodic  Subarrays  with  Random  Displacements 

A  random  array  has  been  suggested  for  grating  lobe  reduction  at  the 
expense  of  higher  side  lobes  in  a  WSN.  We  would  like  to  extend  the  idea  from  using  a 
single  isotropic  antenna  to  a  subarray  at  each  node  location.  The  random  distribution  of 
subarrays  will  follow  a  certain  PDF,  for  example,  uniform,  Gaussian,  truncated  Gaussian, 
Rayleigh,  truncated  Rayleigh,  etc. 

By  applying  the  same  concept  of  pattern  multiplication  from  periodic 
distributed  subarrays,  we  can  easily  find  the  average  radiation  pattern  from  the  single 
subarray  pattern  and  the  randomly  distributed  construction  array  pattern.  The  average 
array  radiation  pattern  formulas  can  be  multiplied  by  a  single  subarray  pattern  (if  all 
subarrays  are  identical),  and  the  final  pattern  for  the  random  distributed  subarrays  can  be 
obtained.  However,  in  order  to  use  the  random  properties  effectively,  the  number  of 
subarrays  has  to  be  large  enough. 
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b.  Periodic  Subarrays  with  Random  Sizes  and  Random 
Displacements 

Subarraying  within  a  single  contiguous  array  has  been  applied  to  minimize 
the  number  of  true  time  delay  units  for  broadband  application.  Ideally,  a  true  time  delay 
unit  at  each  element  is  preferred  for  a  broadband  array.  However,  practical  limitations 
like  cost,  weight  and  complexity  make  it  impossible  to  be  realized.  Therefore, 
subarraying  becomes  the  option  to  improve  the  inherent  bandwidth  limitation  of  a  phased 
array.  The  design  with  contiguous  equal  size  subarrays  results  in  high  side  lobes,  grating 
lobes,  gain  reduction  and  squinting  in  the  pointing  angle  [43].  A  typical  periodic  array 
model  with  random  subarray  sizes  and  subarray  centers  was  presented  in  [47]  and  is 
shown  in  Figure  13. 
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Figure  13.  Generic  array  with  902  elements  and  32  subarrays  and  corresponding 

subarray  centers  (from  [47]). 


For  an  example  of  the  reduction  possible,  see  Figure  14.  A  radiation 
pattern  of  a  linear  array  of  128  elements,  equally  divided  into  32  contiguous  subarrays  of 

o 

four  elements  each,  is  depicted.  For  a  steering  angle  of  40  ,  grating  lobes  can  be  clearly 
seen  due  to  the  periodicity  in  the  array.  The  desired  (average)  radiation  pattern  of  the 
same  array,  but  randomly  divided  and  distributed,  is  shown  in  Figure  15.  The  elimination 
of  the  grating  lobes  and  the  reduction  of  the  side  lobes  are  expected,  since  the  periodicity 
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in  the  array  was  removed.  Also  these  are  applications  where  the  number  of  elements  in 
the  subarrays  is  different  and  could  be  selected  randomly  from  a  discrete  distribution. 


Figure  14.  Pattern  of  a  uniformly  divided  array,  128  elements  divided  into  32 

o 

subarrays  of  four  elements  each,  steering  to  40  (from  [43]). 


Figure  15.  Average  pattern  of  a  randomly  divided  array  (from  [43]). 


32 


D. 


SUMMARY  AND  CONCLUSIONS 


Basic  formulas  for  the  DDSA  radiation  pattern  were  presented  in  this  chapter. 
Both  periodic  and  aperiodic  or  random  arrangements  were  considered,  and  general 
equations  are  derived.  Several  methods  of  potential  interest  for  grating  lobe  suppression, 
including  subarray  rotation,  multiplicative  beamforming  and  randomness,  were  discussed 
in  this  chapter  as  well. 
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III.  PROBLEM  STATEMENT  AND  SOLUTION  APPROACH 


A.  DISTRIBUTED  SUBARRAYS 

The  capabilities  of  the  digital  distributed  subarray  antennas,  such  as  higher 
angular  resolution  and  gain  compared  with  using  single  arrays  individually,  were 
mentioned  in  Chapter  I.  Some  challenges  that  accompany  DDSAs  include  the  grating 
lobe  problem,  calibration,  and  time  and  frequency  synchronization.  Among  the  three,  the 
calibration  (including  error  correction)  and  position  location  problems  are  an  active  area 
of  investigation  (e.g.,  for  distributed  sensor  networks,  distributed  arrays,  etc.)  and  several 
methods  have  been  presented  in  the  literature  to  address  the  problems  [48-52].  Timing 
and  frequency/phase  synchronization  are  required  for  coherent  DDSA  operation  as 
well  [53-56].  For  our  applications  of  interest,  all  arrays  are  on  a  common  structure  (e.g.,  a 
ship  or  a  building),  and  thus  distributing  reference  signals  is  relatively  straightforward. 
The  mitigation  of  grating  lobes,  the  presence  of  which  significantly  affects  the 
performance  of  a  DDSA,  has  received  relatively  little  attention  in  the  past  [4,  57,  58]. 

Grating  lobes  for  a  periodic  DDSA  come  from  the  widely  spaced  subarray 
centers.  The  only  method  to  lower  the  grating  lobes  on  the  transmitting  side  is  to  break  up 
the  periodicity  of  the  subarray  centers.  However,  on  the  receiving  side  there  are  more 
options  available  because  digital  signal  processing  can  be  applied. 

B.  SOLUTION  APPROACHES 

For  our  research,  individual  subarrays  should  be  able  to  operate  independently, 
which  means  that  subarray  grating  lobes  will  not  exist  (i.e.,  the  element  spacing  should  be 
less  than  A/2  for  the  elements  in  every  subarray).  Hence,  thinning  within  the  subarrays 
will  not  be  an  option  for  our  approach.  From  pattern  multiplication,  we  know  that  the 
radiation  pattern  of  a  periodic  distributed  subarray  can  be  expressed  in  terms  of  the 
multiplication  of  a  single  subarray  factor  and  the  construction  factor  (the  array  factor, 
composed  of  the  centers  of  each  subarray).  Grating  lobes  come  from  the  construction 
factor  because  of  the  large  spacing  between  subarrays.  The  objective  is  to  eliminate  or 
suppress  the  grating  lobes  from  the  construction  factor,  so  they  are  at  the  desired  side 
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lobe  level.  From  the  formulas  presented  in  Chapter  II,  we  know  that  the  total  array  factor 
is  the  product  of  the  individual  subarray  pattern  and  the  construction  pattern.  Two  grating 
lobe  reduction  approaches  are  apparent.  One  is  to  reduce  the  construction  lobes  by 
breaking  the  periodicity.  The  other  is  to  suppress  the  grating  lobes  with  the  subarray 
pattern. 

We  first  look  at  some  traditional  methods  for  grating  lobe  suppression  of  a  sparse 
array.  Few  of  the  traditional  methods  have  been  applied  to  distributed  subarrays; 
therefore,  the  methods  are  extended  to  DDSAs,  and  their  effectiveness  is  examined  by 
numerical  simulations.  Due  to  limited  effectiveness  in  suppressing  DDSA  grating  lobes 
by  traditional  methods,  we  investigate  a  new  approach  which  utilizes  the  power  of  digital 
signal  processing  to  eliminate  the  grating  lobes  and  further  lower  the  side  lobes  on  the 
receiving  pattern.  It  is  demonstrated  that  a  combination  of  the  traditional  methods  on  the 
transmitting  side  together  with  the  digital  processing  on  the  receiving  side  provide 
significant  two-way  pattern  improvement. 

1.  Traditional  Solutions 

a.  Sequential  Subarray  Rotation 

As  described  in  Chapter  II,  subarray  rotation  is  capable  of  suppressing  the 
grating  lobes  of  a  sparse  periodic  array.  The  main  idea  is  to  relocate  the  grating  lobes  in 
direction  cosine  space  so  that  the  magnitudes  of  the  lobes  will  not  accumulate.  However, 
subarray  rotation  complicates  the  hardware  design  and  introduces  other  problems,  such  as 
gain  loss  if  the  rotation  angles  are  too  large.  The  grating  lobe  suppression  efficiency 
degrades  as  the  subarray  gap  size  increases  as  well. 

b.  Multiplicative  Beamforming 

Multiplicative  beamforming  on  the  receiving  side  was  introduced  in 
Chapter  II.  The  grating  lobes  can  be  suppressed  by  multiplication  of  the  main  and 
auxiliary  array  outputs.  Careful  design  is  needed  to  place  the  nulls  of  the  auxiliary  array 
on  top  of  the  grating  lobes  of  the  main  array.  Multiplicative  bcamforming  can  only  be 
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used  on  the  receiving  side  and  the  overall  SNR-gain  tends  to  be  lower  due  to  the 
cancellation  of  energy  in  the  grating  lobes. 

c.  Aperiodic  or  Random  Subarray  Sizes 

Random  subarraying  has  been  proposed  for  a  contiguous  subarray  antenna 
to  reduce  the  grating  lobes  due  to  the  subarray  steering  or  when  the  phase  center 
distances  between  subarrays  are  too  wide.  Randomizing  the  subarray  sizes  can  only  lower 
the  grating  lobe  to  a  moderate  level  if  the  subarray  centers  are  still  periodic. 

d.  Aperiodic  or  Random  Displacement 

An  aperiodic  or  random  array  has  no  grating  lobes  because  there  is  no 
strong  periodicity.  For  the  DDSA  case,  aperiodic  or  random  displacement  of  identical 
subarrays  is  an  option  for  lowering  the  grating  lobes.  However,  as  in  the  case  of  aperiodic 
or  random  subarray  sizes,  the  degree  of  grating  lobe  reduction  is  only  moderate. 

2.  Overview  of  Proposed  Solution 

a.  Combination  of  Fundamental  Solutions 

Fundamental  solutions  provide  at  most  moderate  (several  dB)  of  grating 
lobe  suppression.  Multiplicative  beamforming  can  only  be  applied  on  the  receiving  side. 
Subarray  rotation  works  better  for  small  subarray  spacings.  A  combination  of  these 
fundamental  methods  might  provide  better  results  compared  with  their  individual  use. 
Manufacturing  complexity  is  another  important  factor  that  needs  to  be  taken  into  account 
when  designing  a  combined  solution. 

b.  Virtual  Filling  Method 

The  virtual  filling  method  can  only  be  used  on  the  receiving  side.  The  idea 
is  to  fill  the  gaps  of  a  DDSA  with  virtual  elements  so  that  no  grating  lobes  appear  in  the 
visible  region.  After  filling  the  gaps  between  subarrays  virtually,  the  “filled”  array  can  be 
treated  as  contiguous,  and  therefore  some  advantages  that  come  with  a  larger  contiguous 
aperture  array  are  obtained. 
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c.  Combination  of  Fundamental  Solutions  with  Virtual  Filling 
Method 

Since  the  virtual  filling  method  is  applicable  when  receiving  only,  a 
combination  of  fundamental  solutions  is  employed  on  the  transmitting  side.  For  a  radar 
application,  a  two-way  pattern  is  generated  by  multiplying  the  transmitting  pattern  with 
the  receiving  pattern.  Improvement  in  terms  of  lower  side  lobes  and  grating  lobes  is 
possible  after  the  pattern  multiplication. 
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IV.  GRATING  LOBE  SUPPRESSION  WITH  CONVENTIONAL 

SOLUTIONS 


In  this  chapter  we  apply  the  conventional  methods  to  DDSAs  and  examine  them 
individually  to  evaluate  their  effectiveness. 

A.  MULTIPLICATIVE  BEAMFORMING 

The  DDSA  model  shown  in  Figure  16  was  created  for  investigating  multiplicative 
beamforming  on  the  receiving  side.  Arrays  1  through  5  are  4  by  7  rectangular  subarrays 
with  element  spacings  equal  to  half  a  wavelength.  Spacings  between  centers  for  subarrays 
1  through  5  is  six  wavelengths.  Two  auxiliary  subarrays  denoted  6  (8  by  1)  and  7  (1  by 
12)  are  used  to  do  the  pattern  multiplications  for  <f  =  90  and  ^  =  0  planes,  respectively. 
The  center  locations  of  subarrays  6  and  7  are  ( 45  A ,  \0A )  and  ( 60 A ,  20 A ),  respectively. 
Distances  between  the  primary  DDSA  (1  through  5)  and  auxiliary  subarrays  (6  and  7)  are 
not  critical. 
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Figure  16.  Multiplicative  beamforming  DDSA  model. 
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Unlike  traditional  multiplicative  beamforming  using  one  regular  and  one  sparse 
array,  we  would  like  to  have  subarrays  that  can  work  individually  or  collaboratively  if 
needed.  Therefore,  sparse  arrays  are  not  an  option  for  our  model. 

The  radiation  pattern  of  the  primary  DDSA  (arrays  1  through  5)  is  shown  in 
Figure  17  for  the  <p  =  0  plane.  By  multiplying  the  pattern  of  the  primary  DDSA  with  the 
pattern  of  the  auxiliary  subarray  (number  7),  grating  lobes  are  suppressed  dramatically  as 
seen  in  Figure  18. 


Figure  17. 


Primary  DDSA  pattern  (subarrays  1  through  5)  in  the  $  =  0  plane. 
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Figure  18.  Multiplicative  beamforming  pattern  of  arrays  1  through  5  and  auxiliary 
number  7  with  Taylor  amplitude  taper  ( W  =5,  SLL=— 26  dB)  on  arrays  1 
through  5  in  the  (j>  =  0  plane. 


The  same  methodology  can  be  applied  to  the  elevation  pattern  (</>  =  90  ).  The 
original  j>  =  90  pattern  and  the  multiplicative  processed  pattern  are  shown  in  Figure  19 
and  Figure  20,  respectively.  Not  only  has  the  side  lobe  level  been  suppressed,  but  the 
mainbeam  beamwidth  is  reduced. 
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Figure  19.  Primary  DDSA  pattern  of  arrays  1  through  5  in  the  (/  =  90  plane. 
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Figure  20.  Multiplicative  beamformed  pattern  of  arrays  1  through  5  and  auxiliary 
number  6  with  a  Taylor  amplitude  taper  ( W  =5,  SLL=— 30  dB)  on  6  in  the 

(f  =  90  plane. 
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B. 


CONVENTIONAL  ARRAY  GEOMETRY  SOLUTIONS 


For  attacking  the  grating  lobe  problem  on  the  transmitting  side,  introducing 
randomness  is  one  of  the  possible  approaches.  We  consider  subarrays  with  random 
displacements,  rotations  and  sizes.  The  goal  is  to  reduce  the  grating  lobes  that  come  from 
the  periodic  spacing  between  subarrays.  Simulation  models  shown  in  this  section  are  the 
ones  with  the  lowest  grating  lobes  based  on  one  hundred  Monte  Carlo  trials  of  which  the 
maximum  grating  lobes  are  identified  for  all  trials.  It  should  be  noted  that  for  other 
numbers  of  trials,  the  best  arrangement  might  be  different. 

Suppose  we  would  like  to  design  a  periodic  distributed  digital  subarray  antenna. 
Each  subarray  is  a  periodic  array  with  ten  elements  along  x  and  four  elements  along  y, 
with  an  element  spacing  of  ?J 2  in  both  dimensions.  Subarrays  are  separated  by  1.5  A . 
From  Figure  21  and  Figure  22  can  be  seen  the  physical  layout  and  radiation  pattern, 
respectively.  Grating  lobes  occur  due  to  the  periodic  gaps  between  subarrays. 
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Figure  21 .  Physical  layout  of  a  periodic  distributed  linear  array  composed  of  ten 
identical  planar  subarrays  whose  centers  are  equally  spaced.  The  gap  between 

subarrys  is  1.5/1. 
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Figure  22.  Radiation  pattern  of  the  periodic  distributed  subarray  antenna  shown  in 

Figure  21. 

1.  Sequential  Subarray  Rotation 

One  can  use  subarray  rotation  to  mitigate  or  suppress  the  grating  lobes.  In  this 

o 

example,  each  subarray  is  rotated  3  sequentially,  as  shown  in  Figure  23.  From  Figure 
24  one  can  see  the  effect  of  subarray  rotation,  which  is  that  the  grating  lobes  are 
suppressed. 

A  Taylor  amplitude  taper  (JT=  5,  SLL=— 20  dB)  was  used  to  lower  the  side  lobes 
as  shown  in  Figure  25. 
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Figure  23.  Physical  layout  of  the  periodic  distributed  subarray  with  subarray  rotation. 
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Figure  24.  Radiation  pattern  of  the  array  shown  in  Figure  23. 
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Figure  25.  Radiation  pattern  of  the  array  shown  in  Figure  23  with  Taylor  amplitude 

tapering  ( W  =5,  SLL=— 20  dB). 

2.  Random  Subarrav  Sizes 

Next  we  consider  randomly  changing  the  subarray  sizes,  yet  keeping  the  total 
number  of  elements  along  x  and  y  the  same  as  in  the  previous  example  for  comparison 
puiposes.  A  Rayleigh  distribution  is  arbitrarily  chosen  as  the  PDF  for  the  randomization 
process.  However,  it  is  worth  noting  that  the  larger  the  differences  between  the  sizes  of 
the  subarrays  the  better  the  grating  lobe  suppression  capability,  which  means  it  is 
preferable  to  have  subarray  sizes  distributed  over  a  wider  range.  The  physical  layout  is 
shown  in  Figure  26.  By  doing  subarray  size  randomization,  grating  lobes  drop  to  —16  dB 
as  shown  in  Figure  27. 
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Figure  26.  Physical  layout  of  the  random  sized  distributed  subarrays  with  total 
number  of  elements  along  x  and  y  the  same  as  for  the  periodic  case. 
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Figure  27.  Radiation  pattern  of  the  array  shown  in  Figure  26. 


3.  Random  Subarrav  Displacement  and  Sizes 

The  next  step  is  to  add  random  positions  for  each  subarray  on  top  of  their  random 
sizes.  From  Figure  28  and  Figure  29,  one  can  see  the  physical  layout  and  the  radiation 
pattern  of  this  array,  respectively.  As  can  be  seen  in  Figure  29,  the  peak  grating  lobe  level 
has  been  lowered  by  around  1.5  dB  by  doing  random  locations. 
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Figure  28.  Physical  layout  of  the  array  with  random  subarray  sizes  and  random 

subarray  locations. 
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Figure  29.  Radiation  pattern  of  the  array  shown  in  Figure  28. 

4.  Combination  of  Random  Subarray  Displacements,  Sizes  and 
Sequential  Rotations 

As  a  final  step  we  can  use  the  subarray  rotation  idea  from  periodic  distributed 
subarrays.  We  choose  sequential  rotations  instead  of  random  rotations  for  two  reasons. 
First,  sequential  rotations  are  guaranteed  to  move  grating  lobes  to  desired  locations  in 
direction  cosine  space.  Random  rotations  might  move  them  to  nearby  locations  which  can 
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reinforce  the  grating  lobes.  Secondly,  sequential  rotation  is  favorable  in  terms  of 
manufacturing  complexity. 

The  layout  of  the  final  array  with  random  sizes,  position,  and  sequential  rotations 

o 

(2  sequentially)  is  shown  in  Figure  30.  The  radiation  pattern  can  be  seen  in  Figure  31. 
The  peak  grating  lobe  is  approximately  —19  dB,  which  has  only  1  dB  of  improvement 
compared  to  the  non-rotated  case.  Therefore,  this  added  complexity  is  likely  not  worth 
the  effort.  After  applying  the  Taylor  amplitude  tapering  (aF  =  5,  SLL=— 20  dB),  the  final 
pattern  is  shown  in  Figure  32. 


Figure  30.  Physical  layout  of  the  subarray  with  random  subarray  sized  and  random 
subarray  locations  and  sequential  subarray  rotations. 
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Figure  3 1 .  Radiation  pattern  of  the  array  shown  in  Figure  30. 
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Figure  32.  Radiation  pattern  of  the  array  shown  in  Figure  30  with  Taylor  amplitude 
tapering  (W  =5,  SLL=— 20  dB)  in  the  (f  =  0  plane. 

In  Figure  33,  we  summarize  the  effectiveness  of  the  conventional  methods  and 
their  combinations  on  grating  lobe  suppression  for  this  specific  arrangement  without 
amplitude  tapering.  The  improvements  shown  are  computed  relative  to  the  periodic 
DDSA.  For  concise  presentation,  we  define  the  shorthand  notation  of  each  method  as 
below: 

P:  periodic  DDSA. 

RD:  random  subarray  displacement. 

SR:  sequential  subarray  rotation. 

RS:  random  subarray  sizes. 

In  Figure  33,  it  can  be  seen  that  the  combination  of  methods  has  the  greatest 
improvement  relative  to  using  the  conventional  methods  individually.  However,  it  is 
worth  noting  that  an  optimization  process  is  needed  for  each  new  DDSA  arrangement  to 
achieve  maximum  grating  lobe  reduction. 
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Figure  33.  Summary  of  the  effectiveness  on  grating  lobe  suppression  using 
conventional  approaches  individually  and  in  combination. 


C.  SUMMARY  AND  CONCLUSIONS 

In  this  chapter,  several  conventional  methods  for  grating  lobe  suppression  were 
examined.  To  further  improve  the  grating  lobe  suppression  results,  a  combination  of  the 
conventional  methods  was  introduced  and  simulated.  As  can  be  seen  from  the  simulation 
results,  each  method  provides  some  additional  grating  lobe  suppression.  By  combining 
the  methods,  grating  lobes  are  further  lowered  by  approximately  8  dB  (from  -1 1  dB  to 
—19  dB).  However,  disadvantages  such  as  polarization  loss,  hardware  complexity,  gain 
loss  and  limited  suppression  ability  for  large  separations  restrict  the  use  of  these  methods, 
and  therefore,  tradeoffs  need  to  be  made  accordingly. 
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V.  GRATING  LOBE  SUPPRESSION  WITH  VIRTUAL  FILLING 


Complete  digital  control  of  amplitude  and  phase  at  the  element  level  of  an  array 
allows  great  flexibility  in  beamforming.  Modern  radar  and  communications  systems  are 
incorporating  phased  arrays  with  wider  bandwidths,  allowing  for  the  possibility  that 
several  systems  on  the  same  platform  can  share  arrays.  A  system  that  incorporates 
DDSAs  working  cooperatively  as  a  single  array  (thus  forming  an  array  of  subarrays)  can 
potentially  increase  the  output  SNR  and  provide  better  spatial  resolution  compared  with 
using  the  subarrays  individually.  However,  even  if  the  individual  array  patterns  have  no 
grating  lobes,  conventional  beamforming  with  periodic  subarrays  will  have  an  output 
response  with  grating  lobes,  which  is  unacceptable  for  most  applications. 

In  this  chapter,  we  propose  virtual  filling  of  the  gaps  between  the  subarrays  to 
eliminate  the  grating  lobes  on  the  receiving  side  so  that  the  response  of  a  single  large 
contiguous  array  is  synthesized.  Therefore  no  grating  lobes  will  appear  as  long  as  element 
spacing  within  all  subarrays  is  less  than  one  half  of  the  wavelength.  Furthermore, 
amplitude  tapering  can  be  applied  to  the  synthesized  array  to  reduce  interference  and 
clutter.  The  number  of  virtual  elements  that  can  be  used  to  fill  between  or  extend  outside 
the  real  elements  is  restricted  by  the  accuracy  of  the  estimated  parameters. 

Multiple  signals  that  impinge  on  the  array  can  be  either  desired  (e.g.,  radar  target 
return)  or  undesired  (e.g.,  interference  or  clutter).  If  the  subarrays  are  widely  separated, 
then  closely  spaced  grating  lobes  occur,  and  there  will  be  many  angles  where  the 
undesired  signal  has  a  large  response  as  the  main  beam  is  scanned.  We  would  like  an 
output  response  for  the  synthetic  array  that  mimics  the  response  of  a  contiguous  array  so 
that  the  mainbeam  is  in  the  direction  of  the  desired  signal  and  the  interference  is  in  a  low 
side  lobe.  This  response  can  be  synthesized  by  filling  in  the  gaps  between  the  arrays  with 
virtual  elements.  Here,  methods  are  used  to  estimate  the  virtual  element  weights  from  the 
in-phase  and  quadrature  baseband  signals  received  by  the  real  elements.  The  information 
needed  to  reconstruct  the  contiguous  array  response  is  the  signal  DOAs,  magnitudes  and 
phases. 
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In  Section  A,  the  DDSA  model  is  introduced.  Because  DOA  estimation  is  crucial 
to  synthesizing  the  virtual  element  weights,  Section  B  discusses  how  some  single 
snapshot  (SS)  DOA  algorithms  perform  with  regard  to  this  problem.  In  particular,  the 
modified  MP  algorithm  is  found  to  have  improved  resolution  for  closely  spaced  signals 
relative  to  the  Fast  Fourier  Transform  (FFT)  method.  A  multiple  snapshot  (MS) 
algorithm  is  also  described  in  Section  B.  In  Section  C,  we  present  formulas  for  estimating 
signal  amplitudes  and  phases  based  on  the  DOAs.  The  dependence  on  element  level  SNR 
of  the  synthesized  array  response  is  formulated.  Noise  and  other  errors  are  considered  in 
Section  D. 

A.  DDSA  MODEL 

For  simplicity  we  consider  a  linear  array  of  Ns  identical  subarrays.  Each  subarray 
contains  Nx  elements  that  are  equally  spaced  dx  along  the  x  axis  as  shown  in  Figure  34. 
The  gap  between  subarrays  is  D  =  hdx  where  h  is  an  integer  greater  than  zero  (i.e.,  the 
gap  is  an  integer  multiple  of  the  element  spacing  as  shown  in  Figure  34). 


2 


Figure  34.  Linear  distributed  digital  subarray  model.  Black  filled  dots  are  real 
elements  and  non-filled  dots  are  virtual  elements. 
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If  the  entire  array  is  centered  at  the  origin,  the  location  of  element  p  in  subarray  m 


is 


(34) 

=  P{m,p)dx  (m  =  1,2, ...,  Ns;  p  =  1,2, Nx). 

If  there  are  K  <  NsNx  signals  incident  on  the  array  from  angles  0r 
(r  =  l,2, K)  with  complex  voltages  Vre,a’  the  element  outputs  can  be  expressed  in 
phasor  form  as: 

K 

A(m,  p)  =  I  V  exp (~jkP(m,  p)dx  sin  9  +  jar ) 

r=\  (35) 

=  I(m,p)  +  jQ(m,p). 


x(m,p)  = 


P~ 


N,NX  +(Ns-\)h  +  \ 


+  {m-  \)(NX  +  h) 


B.  DOA  ESTIMATION 

Numerous  DOA  estimation  algorithms  are  available,  but  the  MP  method  performs 
particularly  well  for  single  snapshot  noisy  data.  It  utilizes  singular  value  decomposition 
(SVD)  to  divide  the  matrix  space  into  signal  and  noise  subspaces.  By  discarding  the 
eigenvector  corresponding  to  the  noise  signal,  the  noise  effect  can  be  reduced,  and  hence, 
the  estimation  accuracy  can  be  improved. 

From  measurement  of  the  I  and  Q  at  the  elements,  the  signal  parameters  0r  and 
VreJa'  can  be  estimated.  Thermal  noise  is  accounted  for  by  adding  a  complex  noise  to  the 
A{m,p)  in  Eq.  (35).  The  noise  leads  to  an  error  in  the  parameter  estimates,  which  in  turn 
results  in  a  distortion  of  the  synthesized  antenna  response. 

1.  Matrix  Pencil  Method 

General  formulas  of  the  classical  MP  method  used  to  estimate  the  DOAs  are 
presented.  A  single  one-dimensional  (1-D)  uniform  linear  array  model  is  shown  in  Figure 
35. 
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Figure  35.  Linear  array  model. 


For  a  linear  array  centered  at  the  origin,  locations  of  each  element  can  be 
represented  similar  to  Eq.  (34)  as: 


(36) 


The  sum  of  K  complex  signals  received  at  each  element  can  be  modeled  as  follows: 

K 


A(P)  =  XK  exp (~jkP(p)dx  sin  0,  +  jar) 


(37) 


=  Kp)+jQ{p) 


where  Vre’a'  is  the  complex  magnitude  associated  with  signal  r ,  0r  is  the  dircction-of  - 
arrival  of  signal  r  with  respect  to  the  z  direction. 

In  the  noisy  environment,  we  use  SVD  directly  with  received  complex  signals  in 
the  form  of  Hankel  structure  matrix  of  dimension  (Nx  —L  +  l)xL,  in  which  L  is  the 
pencil  parameter  chosen  between  N  J2  and  7Vr/3  for  optimum  noise  removal  [23],  The 
Hankel  structure  matrix  A  is 


/ 


^d) 

A(  2) 


A(  2) 
A(  3) 


A(L) 
A(L  + 1) 


A  = 


(38) 


A(NX-L  + 1)  A(NX-L  +  2)  •••  A(NX) 


(Nt-L+\)xL 
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The  SVD  of  A  can  be  represented  as: 

A  =  WTU"  (39) 

where  W  and  U  are  unitary  matrices  whose  columns  are  eigenvectors  of  AA/y  and  A/yA, 
respectively.  T  is  a  diagonal  matrix  containing  the  singular  values  of  A.  If  the  number  of 
signals  (K)  is  known,  we  can  select  the  largest  K  singular  values  and  their  corresponding 
eigenvectors,  and  form  a  new  signal  matrix  W  which  contains  only  the  signal  subspace 
eigenvectors  and  has  most  of  the  noise  effects  removed. 


Next,  by  solving  the  following  eigenvalue  problem,  the  DOAs  can  be  estimated: 


(W,tW2)v  =  Arv  r  =  1, ...,  K 

Wj  =  (Wl//W1)",Wl" 


(40) 


where  v  is  the  eigenvector,  and  A.  is  the  eigenvalue  of  W,+W2 ,  W,  and  W2  are  W  with 
the  last  and  first  row  deleted,  respectively.  The  symbol  f  is  the  Moore-Pcnrose  pseudo 
inverse  as  defined  in  Eq.  (40).  We  then  can  extract  the  DOAs  from  Ar  using  the 
following  equation: 


i  f  Im[ln(A,,)]] 

I  w  y 


(41) 


2.  Modified  Matrix  Pencil  Method  for  DDSA 
a.  Single-Snapshot  MP  Method 

We  propose  an  extension  of  the  MP  method  that  is  tailored  to  the  DDSA 
by  arranging  the  Hankel  matrices  of  each  subarray  from  top  to  bottom  sequentially.  Let 
Ym  be  the  Hankel  matrix  for  subarray  tn 


Y 


m 


A(m,  1) 

A(m,Nx-L  + 1) 


A(m,L) 

A(m9Nx) 


\ 


/[{Nx-L+\)\xL 


(42) 
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so  that  for  a  single  snapshot 


lss 


'  '[{Nx-L+\)Nt]*L 


(43) 


b.  Multiple  Single-Snapshot  MP  Method 

Due  to  rapidly  advancing  technology,  receivers  are  capable  of  high 
sampling  rates.  Therefore,  a  multiple  single-snapshot  (MSS)  MP  method  is  relatively 
straightforward,  and  it  is  shown  in  this  section  to  have  good  thermal  noise  removal 
capability  in  terms  of  DOA  estimations.  Based  on  the  simulation  results,  thirty  single 
snapshots  will  be  enough  for  acceptable  DOA  estimation.  Using  more  single  snapshots 
improves  the  estimation  accuracy  at  the  expense  of  longer  waiting  time  and  a  requirement 
for  more  computational  power. 

A  multiple  single-snapshot  method  is  based  on  the  single-snapshot 
formulas  described  in  the  previous  section.  By  averaging  the  estimated  DOAs  from 
snapshots,  the  thermal  noise  effects  on  DOA  estimations  can  be  reduced  dramatically. 


c.  Multiple  Snapshots 

In  principle,  MP  requires  only  a  single  snapshot,  but  it  can  be  extended  to 
multiple  snapshots,  thus  resulting  in  a  lower  root  mean  square  error  (RMSE)  [28]. 


The  multiple  snapshot  MP  can  be  considered  as  a  concatenation  of 
multiple  columns  of  single  snapshot  MP.  If  Ym  h  is  the  Hankel  matrix  for  snapshot  b 

(b  =  1,2,.  ..,5)  of  subarray  tn 


Y  ,  = 

ttlyO 


A„(m,  1)  ...  Ah{m,L) 


A„(m,Nx-L  + 1)  Ah{m,Nx)){N'_L^L 


(44) 


58 


then  for  multiple  snapshots 


y  = 

1  A/S 


.V(  ,1 


\ 


S[(Nx-L+l)Na]x(LB) 


(45) 


Note  that  computational  time  increases  with  the  number  of  snapshots. 
Some  advantages  of  using  multiple  snapshots  are  to  stabilize  the  DOA  estimation  for  a 
small  number  (one  or  two)  of  subarrays  cases  at  lower  SNR  (0  to  4  dB).  The  RMSE  of 
DOA  estimation  tends  to  decrease  when  the  number  of  snapshots  increases.  However,  in 
terms  of  DOA  estimation,  the  average  of  multiple  single-snapshots  will  provide  more 
accurate  results  compared  to  multiple  snapshot  case  as  in  Eq.  (45).  This  is  because  the 
noise  effects  will  be  reduced  again  by  averaging. 


d.  Simulation  Results 

After  the  matrices  for  the  DDSA  are  formed,  the  standard  MP  procedure 
for  finding  the  DO  As  of  signals  in  the  noisy  environment  is  applied  [32]. 

The  improved  performance  of  the  angle  estimates  from  the  modified  MP 
method  was  verified  using  a  Monte  Carlo  simulation  with  100  trials.  First,  a  signal  is 
incident  from  30  with  phase  7r/5  radian  onto  a  DDSA  comprised  of  five  eight-element 
linear  arrays  with  an  element  spacing  of  0.42/1 .  The  spacing  between  subarrays  is  3.36/1 . 
The  advantage  of  using  the  modified  MP  method  can  be  observed  in  Figure  36.  It  can  be 
seen  that  for  one  subarray  (i.e.,  small  NSNX )  at  low  SNR,  the  RMSE  of  the  DOA  is  high. 
As  expected,  employing  more  subarrays  makes  the  estimates  much  more  accurate. 
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Figure  36.  RMSE  of  DOA  for  one  signal  vs.  SNR  per  element  for  various  numbers  of 
subarrays.  Each  subarray  has  eight  elements  spaced  0.42/1 ,  and  gaps  are  equal 

to  the  subarray  size. 


Next,  we  increase  the  number  of  elements  in  the  DDSAs  from  eight  to 
eighty.  The  gaps  are  now  80^/  =  33.6A  long.  Also  the  number  of  signals  is  increased  to 
two:  one  from  30  and  a  second  from  31  .  The  phases  of  the  signals  are  nt 5  and  -4/r/5, 
respectively.  From  the  curves  in  Figure  37  and  Figure  38,  it  can  be  seen  that  increasing 
the  number  of  DDSAs  provides  much  smaller  RMSE  even  at  low  SNR  (0  dB  )  for  the  two 
signals.  If  the  FFT  method  [59]  were  used  at  the  subarray  level,  the  beamwidth  would  be 
too  wide  to  resolve  the  two  closely  spaced  signals.  As  for  the  single  target  case,  using 
more  subarrays  in  the  processing  results  in  much  more  stable  and  accurate  results,  which 
will  be  crucial  for  effective  “filling.” 
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Figure  37.  RMSE  of  DOA  for  signal  1  vs.  SNR  per  element  for  various  numbers  of 
subarrays.  Each  subarray  has  eighty  elements  spaced  0.42/1 ,  and  gaps  are  equal 

to  the  subarray  size. 
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Figure  38.  RMSE  of  DOA  for  signal  2  vs.  SNR  per  element  for  various  numbers  of 
subarrays.  Each  subarray  has  eighty  elements  spaced  0.42/1 ,  and  gaps  are  equal 

to  the  subarray  size. 
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C.  VIRTUAL  FILLING  METHOD 
1.  Introduction 

The  sum  of  K  complex  signals  from  the  source  Vrej°r  from  angle  6y  at  the  output 
of  the  elements  of  subarray  m  can  be  written  as 

K 

A(m,p)  =  exp(- yXPO,  p)dx  sin  dr)txp(jar)  (p  =  \,...,NX )  (46) 

/=1 

Casting  these  in  matrix  form  gives 

K  V  =  Am  (47) 

where 

-  jkP(m,l  )d  sin  Qx  -  jkP(m,\  )Jsin  9K 

t  •  •  •  L 

-  jkP(m,Nx  )</sin^,  -  jkP{m,N x  )ds\nOK 


Note  that  adding  subarrays  increases  the  total  number  of  elements  and  thus  the  number  of 
signals  that  can  be  handled. 
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Now  the  estimated  signal  magnitudes  V  and  phases  ar  can  be  used  to  create 
virtual  complex  data  to  “fill”  the  gaps  between  subarrays.  The  location  of  virtual  element 
g  in  gap  /  (between  subarrays  i  and  /+1)  is 


x(i,g)=  g- 


(Ns  -  2)Nx  +  (Ns  - \)h  +  1 


+  (i-l)(Nx+h) 


d. 


(52) 


=  Z(i,g)dx  ( g  =  1,2,..., A;  i  =  l,2,...,Ns  - 1) 


The  complex  data  for  filling  is  given  by  the  same  formula  as  for  the  real  data  in  Eq.  (46) 

B(i,g)  =  I  Vrexp (~jkZ(i,g)dx  sin^.)exp(y«,.)-  (53) 

r= 1 


Combining  the  real  and  virtual  data  gives  the  response  of  the  synthesized  array 


F  = 


(■ A|  W  (Bi  )ix* ' '  ‘  (B.y,-i  )lx/,  )lxlV 


(54) 


'  J\x[NsNx+(Ns-l)h) 


where  the  B  partitions  are  comprised  of  the  tenns  given  by  Eq.  (53).  F  represents  the 
complex  outputs  of  the  synthesized  array  in  the  K  signal  directions.  Multiplying  by  the 
desired  beamforming  weights  and  summing  gives  the  array  response. 


2.  Simulation  Results 

Consider  a  five-subarray  DDSA  with  30  elements  in  each  subarray  and  an 
element  spacing  of  0.42/1.  The  subarray  length  is  12.6A.  The  gaps  are  also  (arbitrarily) 
set  to  12.6/1.  One  unit  amplitude  signal  is  incident  from  0  with  a  phase  of  7r/5.  A 
second  interference  signal  is  coming  in  at  2.3  with  a  phase  -4/r/5.  The  pattern  is  shown 
in  Figure  39  as  the  weights  are  changed  to  scan  the  main  beam  in  a  region  of  direction 
cosine  space  (sin#) .  A  20  dB  Taylor  amplitude  distribution  is  applied.  As  can  be  seen, 
the  high  response  of  the  interfering  signal  that  occurs  at  grating  lobe  locations  has  been 
eliminated.  The  synthesized  array  response  in  the  direction  of  both  signals  is  the  same  as 
that  of  a  contiguous  array. 
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Figure  39.  Comparison  of  original  and  synthesized  antenna  responses  of  5  subarrays 

o  o 

each  with  30  elements,  for  signals  (noiseless)  incident  from  0  and  2.3  . 

The  average  synthesized  array  response  for  the  same  case  shown  in  Figure  39  can 
be  seen  in  Figure  40,  but  with  an  SNR  per  element  of  6  dB  (single  snapshot).  The  results 
of  a  Monte  Carlo  simulation  of  thirty  trials  (equivalent  to  thirty  single  snapshots)  using 
Eq.  (43)  for  the  five-subarray  two-signal  case  used  to  generate  Figure  40  are  summarized 
in  Table  1. 


64 


Figure  40.  Comparison  of  original  and  synthesized  antenna  response  after  virtual 
filling  for  an  element  level  SNR  of  6  dB.  Taylor  amplitude  taper  (W  =5,  SLL= 

—20  dB)  has  been  applied. 


Table  1.  Exact  and  estimated  signal  parameters  for  virtual  processing. 


Parameter 

Actual 

Estimated 

(mean) 

Estimated 

(variance) 

Signal  1 

Signal  2 

Signal  1 

Signal  2 

Signal  1 

Signal  2 

DOA 

(degrees) 

0 

2.3 

0.0356 

2.2965 

0.1808 

0.2878 

Amplitude 

(Volts) 

1 

1 

0.9630 

0.9533 

0.002 

0.0034 

Phase 

(radians) 

0.6283 

2.5133 

0.6371 

2.5045 

0.0018 

0.002 

The  effectiveness  of  the  virtual  approach  relies  on  having  accurate  angle  estimates 
for  the  signals  of  interest,  which  in  turn  requires  a  high  effective  array  SNR  (i.e.,  large 
N  N x  if  the  SNR  per  element  is  low;  high  SNR  if  the  total  number  of  elements  NSNX  is 
small).  The  SNR  per  element  can  be  increased  by  adding  a  low  noise  amplifier  or 
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increasing  the  gain  of  the  array  element.  It  is  also  possible  to  improve  the  DOA  estimates 
with  multiple  snapshots  or  by  averaging  multiple  single  snapshots. 

D.  NOISE  AND  OTHER  ERRORS 

For  a  real  application,  factors  such  as  noise  and  errors  cause  performance 
degradation  of  a  system.  Thermal  noise  appears  in  almost  all  receivers  and  is  commonly 
considered  as  the  primary  error  source  in  the  receiver.  Low  noise  amplifiers  are  usually 
used  to  increase  the  SNR.  Other  errors  including  array  element  positioning  error,  mutual 
coupling  error,  manufacturing  imperfections  and  calibration  errors  are  considered  as 
fixed  errors  that  can  be  measured  in  advance  and  possibly  be  compensated.  We  model  the 
thermal  noise  at  the  element  level  by  specifying  the  element  SNR.  Not  only  is  it  random 
between  elements,  but  it  is  also  random  from  snapshot  to  snapshot.  On  the  other  hand, 
fixed  errors  are  modeled  as  random  variables  that  are  fixed  between  snapshots/trials. 

Next,  a  five-subarray  DDSA  model  is  used  to  examine  the  effects  of  fixed  errors 
to  the  DO  As  estimation  and  filling  method.  Each  subarray  is  comprised  of  ten  elements 
with  element  spacings  equal  to  0.45 A.  Subarray  center  distances  are  \0A.  Fixed  errors 

o  o  o 

are  uniformly  distributed  from  -21  to  21  (root  mean  square  (RMS)  values  from  0  to 

o 

12.1  ),  and  the  SNR  is  varied  from  6  dB  to  21  dB  at  each  element.  Two  signals  with 
equal  magnitude  and  non-coherent  phases  (  tt/5  and  -4k  15  )  from  DO  As  of -10  and 

o 

15  relative  to  broadside  are  impinging  on  the  DDSA. 

The  receiving  pattern  for  6  dB  SNR  for  the  ideal  contiguous  array  of  the  same 
aperture  size  as  DDSA,  original  DDSA  and  virtual  filling  method  are  compared  in  Figure 
41.  A  huge  improvement  in  terms  of  grating  lobes  and  sidelobe  suppression  can  be 
observed  after  applying  the  virtual  filling  method. 
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Main  Beam  Scan  Angle  in  sin(0) 


Figure  41.  Pattern  comparison  of  contiguous,  original  DDSA  and  virtual  filling. 
Assuming  no  fixed  errors  and  with  6  dB  SNR  at  each  element.  Taylor 
amplitude  taper  (7T=  5,  SLL=— 30  dB)  has  been  applied. 


In  order  to  quantify  the  effect  of  fixed  errors  on  the  DOA  estimations,  a  plot  that 
compares  the  RMSE  of  the  DOA  versus  RMS  phase  error  for  different  SNR  levels  is 
shown  in  Figure  42.  The  formula  used  to  calculate  the  average  RMSE  of  the  DOA  for  the 
K  signals  is 


RMSE  = 


(55) 


where  9r  is  the  estimated  angle,  and  Gr  is  the  true  angle  of  signal  /*.  Two  conclusions  can 
be  made  based  on  the  information  provided  by  the  plot: 

1.  The  most  efficient  way  to  improve  the  RMSE  of  the  DOA  is  to  increase 
the  SNR  at  each  element. 


o 

2.  For  RMS  phase  errors  less  than  6  which  correspond  to  a  fixed  error  of 

o  o 

-10.5  to  10.5  ,  the  effect  on  DOA  estimation  can  be  ignored  for  all 
SNRs. 
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Figure  42.  RMSE  of  DOA  versus  RMS  phase  errors  from  0  to  12.1  for  different 

element  SNR. 

Many  of  the  sources  of  the  fixed  errors  can  be  compensated  for  by  pre-calculation 
or  pre-measurement,  but  there  will  still  be  some  residual  errors  after  correction.  We 

o 

consider  phase  errors  up  to  21  and  examine  how  they  degrade  the  radiation  pattern. 

o 

Figure  43  has  a  plot  of  the  pattern  of  the  worst  case  (21  fixed  error)  at  6  dB  SNR.  By 
comparing  Figure  43  with  Figure  41,  we  see  that  the  effect  of  the  fixed  error  on  the 
receiving  pattern  is  to  increase  the  side  lobe  level  and  lower  the  main  beam  by  0.6  dB  due 
to  the  increase  of  RMSE  of  the  DOA. 
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Figure  43.  Pattern  comparison  of  contiguous,  original  DDSA  and  virtual  filled  DDSA 

o 

for  21  fixed  error  and  with  6  dB  SNR  at  each  element.  Taylor  amplitude  taper 
(  W  =  5,  SLL=— 30  dB)  has  been  applied. 

We  also  consider  the  effect  of  signal  angles  on  the  RMSE  of  DOA  with  fixed 

o  ° 

errors  uniformly  distributed  between  -21  and  21  and  6  dB  SNR  per  element.  As  can  be 
seen  in  Figure  44,  because  the  projected  aperture  area  decreases  at  large  signal  angles,  the 
RMSE  of  the  DOA  is  increased. 
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Figure  44.  RMSE  of  DOA  versus  RMS  phase  errors  from  0  to  12.1  and  6  dB  SNR 

per  element  for  different  signal  angles. 


E.  SUMMARY  AND  CONCLUSIONS 

Filling  the  gaps  between  arrays  with  virtual  elements  for  the  purpose  of  receiving 
processing  allows  a  synthesized  antenna  response  that  duplicates  a  contiguous  array  in  a 
number  of  directions  that  is  limited  by  the  total  number  of  elements  used  in  the 
processing.  The  synthesized  response  has  no  grating  lobes  and  can  have  low  side  lobes. 
As  a  first  step,  the  signal  amplitudes,  phases  and  DOAs  must  be  extracted  from  the 
element  /  and  Q  samples.  This  data  is  used  generate  I  and  Q  samples  that  would  be 
provided  by  virtual  elements  filling  the  gaps  between  distributed  subarrays.  Low  side 
lobes  and  interference  rejection  were  demonstrated  for  the  virtual  processed  DDSA. 

The  MP  method  was  found  be  to  be  well  suited  for  this  application.  The  MP 
technique  was  extended  to  handle  multiple  subarrays,  for  either  single  or  multiple 
snapshots.  Multiple  snapshots  provide  improved  stability  in  a  low  SNR  situation.  This 
method  can  deal  with  both  coherent  and  non-coherent  signals  and  requires  fewer 
snapshots  for  accurate  DOA  estimation. 
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The  proposed  filling  method  allows  suppression  of  undesired  signals  that  would 
normally  occur  at  grating  lobe  angles.  It  was  shown  that  for  high  SNR  the  array  response 
approaches  that  of  a  contiguous  array  of  the  same  extent. 
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VI.  APPLICATIONS 


A.  THREE-STEP  MODIFIED  MATRIX  PENCIL  METHOD  FOR  DDSA 

From  the  simulation  results  in  Chapter  V,  it  was  shown  that  the  modified  MP  for 
DDSA  provides  improved  performance  compared  to  the  individual  use  of  the  arrays.  The 
overall  estimation  accuracy  can  be  further  improved  by  adding  a  step,  resulting  in  a 
Three-Step  Modified  Matrix  Pencil  method  described  in  this  chapter. 

It  is  well  known  that  the  angular  resolution  of  an  array  antenna  is  inversely 
proportional  to  its  aperture  size.  It  is  true  for  the  MP  method  that  higher  resolution 
requires  larger  array  aperture,  i.e.,  adding  more  radiation  elements.  However,  a  basic 
requirement  for  the  MP  method  is  to  have  a  linear  phase  relationship  between  all  array 
elements.  For  the  DDSA  case,  the  first  workaround  for  the  linear  phase  limitation  is  to 
cascade  the  subarray  Hankel  matrices  as  proposed  in  Chapter  V.  A  second  method 
proposed  here,  is  to  utilize  the  advantage  of  a  large  aperture  and  satisfy  the  linear  phase 
requirement  of  MP  method  to  gradually  improve  the  DOAs  estimation  accuracy.  The 
three  steps  are  delineated  in  the  following  sections. 

1.  Preliminary  DOA  Estimation 

A  preliminary  multiple  single-snapshot  DOA  estimation  using  the  modified  MP 
method  is  the  first  step  of  the  method.  A  fraction  of  the  total  snapshots  is  responsible  for 
this  preliminary  estimation.  For  example,  twenty  single  snapshots,  out  of  a  total  of  forty 
collected,  can  be  used  for  the  preliminary  DOA  estimation. 

2.  Signal  Extraction  and  Virtual  Filling 

Using  the  DOAs  estimated  in  step  one,  we  can  extract  the  signal  parameters  and 
then  virtually  fill  the  gaps  between  subarrays  using  the  estimated  parameters  in  the 
equations  in  Chapter  V,  Section  C. 
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3.  DOA  Estimation  Using  Virtually  Filled  DDSA 

The  next  step  is  to  use  the  “filled”  large  array  for  the  second  improved  DOA 
estimation.  The  filled  array  is  analogous  to  a  large  continuous  array,  and  the  estimation 
results  should  be  close  to  a  real  array  with  the  same  aperture  size  if  the  filling  parameters 
are  accurate  enough.  A  large  aperture  provides  better  angular  resolution.  Therefore  the 
system  performance  in  terms  of  RMSE  should  be  improved  relative  to  the  preliminary 
DOA  estimation. 

4.  Final  DOA  Estimation 

The  final  step  is  to  use  the  improved  estimated  parameters  from  the  previous  step 
to  fill  the  gaps  again.  Since  more  independent  data  points  are  used,  the  average  SNR  at 
each  element  improves.  Higher  SNR  at  each  element  will  contribute  to  a  lower  RMSE  of 
DOA  estimation. 

5.  Simulation  Results 

Two  equal  amplitude  signals  with  phases  7r/5  and  -An! 5  are  impinging  on  a  five- 

o  o 

subarray  DDSA  from  the  angles  of  -20  and  0  relative  to  broadside.  The  spacing 
between  subarray  centers  is  18,4.  Each  subarray  has  20  elements  with  0.45/4  element 
spacing.  The  performance  comparison  of  modified  MP  method  and  three-step  modified 
method,  in  terms  of  RMSE  of  DOA,  is  shown  in  Figure  45.  Due  to  the  increase  of 
aperture  size  of  the  filled  array  1,  the  RMSEs  of  DOAs  are  much  smaller  relative  to  the 
five  subarray  case. 
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Figure  45.  RMSE  of  DOA  for  five  subarrays  (modified  MP)  and  filled  array  1  (three- 

step  modified  method). 


Now  we  use  the  results  (DOAs,  amplitudes  and  phases)  from  filled  array  1  to 
virtually  refill  the  gaps  between  subarrays  and  re-estimate  the  DOA  using  the  MP  method 
for  single  array.  Improvement  in  terms  of  RMSE  can  be  observed  as  seen  in  Figure  46. 
This  is  because  of  the  improved  accuracy  of  the  filling  data  from  the  filled  array  1  that 
lower  the  RMSE  of  DOA  of  filled  array  2. 
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Figure  46.  RMSE  of  DOA  for  filled  array  1  and  filled  array  2. 

6.  Three-Step  Modified  MP  for  Close  Targets 

It  is  known  that  a  larger  aperture  is  the  key  solution  to  resolve  two  closely  spaced 
targets.  We  would  like  to  demonstrate  the  improved  resolution  for  two  close  targets  using 
the  Three-Step  Modified  MP  method  and  compare  the  result  with  that  using  modified  MP 
method.  Two  equal  amplitude  signals  with  phases  n! 5  and  -A7r/5  are  impinging  on  a 

o  ° 

five-subarray  DDSA  from  the  angles  of  -10  and  -4  relative  to  the  broadside.  The 
subarray  spacing  between  centers  is  1 8A  .  Each  subarray  has  20  elements  with  0.45>i 
element  spacing  and  0  dB  SNR  per  element.  A  total  of  40  single-snapshots  are  used  for 
both  modified  MP  method  and  the  Three-Step  Modified  MP  method.  The  radiation  power 
pattern  using  modified  MP  method  is  shown  in  Figure  47.  In  Figure  48,  it  can  be  seen 
that  the  grating  lobes  are  suppressed  further  by  using  the  Three-Step  Modified  MP 
method  as  the  result  of  the  higher  accuracy  of  the  parameters  estimations. 
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Figure  47.  Normalized  power  pattern  of  two  targets  from  -10  and  -4  using  the 

Modified  MP  method.  Taylor  amplitude  taper  ( W  =5,  SLL=— 30  dB)  has  been 

applied. 


°  A° 

Figure  48.  Normalized  power  pattern  of  two  targets  from  -10  and  -4  using  the 
Three-Step  Modified  MP  method.  Taylor  amplitude  taper  (  TT  =  5,  SLL=— 30 

dB)  has  been  applied. 


77 


B.  COMBINATION  OF  RANDOM  SUBARRAY  SIZES  AND  FILLING 

METHOD 

1.  Introduction 

Wc  proposed  the  virtual  filling  method  on  the  receiving  side  to  mimic  a 
contiguous  array  pattern  and  therefore  side  lobe  taper  can  be  applied  to  lower  the  side 
lobes.  However,  the  virtual  filling  method  cannot  be  used  on  the  transmitting  side,  and 
the  grating  lobes  from  a  DDSA  on  the  transmitting  side  can  only  be  addressed  with 
traditional  methods  such  as  random  subarray  displacements,  random/sequential  subarray 
rotations  and  random  subarray  sizes,  as  discussed  in  Chapter  IV.  Wc  choose  random 
subarray  sizes  with  fixed  subarray  gaps  for  our  transmitting  model  because  it  has  the 
following  advantages: 

•  Since  the  gaps  between  subarrays  are  fixed,  it  is  much  easier  to  implement 
compared  to  random  subarray  gaps. 

•  Without  subarray  rotation,  the  virtual  filling  method  can  be  applied  much 
more  easily. 

•  With  random  subarray  sizes  and  fixed  subarray  gaps,  the  subarray  phase 
centers  are  actually  randomized  and  grating  lobes  are  partially  suppressed. 

2.  Simulation  Results 

A  random  sized  20-subarray  DDSA  with  number  of  elements  of  6,  8,  12,  10,  14, 
10,  8,  14,  18,  16,  33,  21,  31,  47,  53,  39,  65,  21,  45,  29  was  first  chosen  on  the  basis  of  the 
lowest  GL  level  from  a  100  Monte  Carlo  simulation  trials  that  randomly  assigned  the 
number  of  elements  of  each  subarray  under  the  constraint  that  the  total  number  of 
elements  in  the  DDSA  was  500.  Element  spacings  are  0.48A ,  and  spacings  between 

-  o  o 

subarrays  are  4.8/1.  Two  signals  of  equal  power  coming  in  from  10  and  -25  relative 
to  broadside  are  used  for  this  simulation.  The  transmitting  pattern  is  noiseless  with  a 
Taylor  amplitude  taper  with  parameters  n  =5,  SLL=— 20  dB.  On  the  receiving  side,  a 
SNR  of  6  dB  at  each  element  and  a  Taylor  amplitude  taper  of  —35  dB  are  used. 
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The  radiation  pattern  of  the  periodic  20-subarray  DDSA  (25  elements  in  each 
subarray)  is  shown  in  Figure  49.  The  high  grating  lobes  due  to  the  large  periodic  spacing 
between  subarrays  are  clearly  seen.  The  pattern  of  the  random  20-subarray  case  is  shown 
in  Figure  50.  The  randomness  has  lowered  the  grating  lobes  but  also  increased  the 
average  side  lobe  level.  We  can  further  lower  the  side  lobe  level  on  the  receiving  side  by 
filling  the  gaps  between  subarrays  so  that  an  amplitude  taper  (n~  =  5,  SLL=— 35  dB)  can 
be  applied  as  shown  in  Figure  51. 


o  o 

Figure  49.  Periodic  DDSA  transmitting  pattern  of  signals  from  10  and  -25  related 

to  the  broadside. 
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Figure  50.  Random  subarray  sizes  DDSA  transmitting  pattern  of  signals  from  10 

o 

and  -25  related  to  the  broadside. 


Figure  51.  Random  subarray  sizes  DDSA  virtual  filling  receiving  pattern  of  signals 

□  o 

from  10  and  -25  relative  to  the  broadside.  A  35  dB  Taylor  distribution  was 

applied  (/f  =  5,  SLL=— 35  dB). 
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c. 


TWO-WAY  PATTERN 


1.  Introduction 

By  applying  pattern  multiplication,  the  normalized  two-way  pattern  is  defined  as 

<*’*)  =  F*ormr,  (#,*)  ■ X  F„armJ0,  </>)  (56) 

where  Fnormj  (#,^)  is  the  normalized  pattern  of  the  transmitting  DDSA  array,  and 
FnorniR  (0,(/)  is  the  normalized  pattern  of  receiving  DDSA  array.  Using  the  two-way 

pattern  multiplication  approach,  we  see  that  the  remaining  DDSA  transmitting  grating 
lobes  can  be  reduced  by  the  pattern  of  the  DDSA  receiving  array. 

The  two-way  beam  pattern  design  of  a  DDSA  [5]  can  lead  to  an  increase  in  the 
hardware  complexity  compared  to  separate  transmitting  and  receiving  systems  and  needs 
to  be  carefully  designed  to  match  the  grating  lobe  and  null  locations  on  the  two  sides. 
Perfect  suppression  of  the  grating  lobes  by  the  nulls  will  only  work  for  a  limited  number 
of  cases.  The  method  we  propose  uses  a  combination  of  random  subarray  sizes  for  the 
transmitting  DDSA  and  virtual  filling  for  the  receiving  DDSA,  as  seen  in  Figure  52; 
therefore,  the  hardware  complexity  is  reduced.  Because  of  the  same  aperture  sizes,  both 
transmitting  and  receiving  patterns  possess  the  same  mainbeam  beamwidth.  The  low 
sidelobe  attribute  of  the  receiving  pattern  can  significantly  reduce  the  grating  lobes  of  the 
transmitting  pattern  after  the  multiplication  has  been  done.  Therefore  no  grating  lobe  and 
null  matching  is  required. 


Real  elements  Virtual  elements 


Figure  52.  The  transmitting  DDSA  model  with  real  elements  (in  blue)  and  the 
receiving  DDSA  model  with  real  and  virtual  elements  (in  blue  and  red 

respectively). 
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2. 


Simulation  Results 


For  comparison  purposes,  a  two-way  pattern  of  the  periodic  DDSA  is  shown  in 
Figure  53.  The  two-way  pattern  shown  in  Figure  54  is  generated  by  multiplying  the 
transmitting  pattern  shown  in  Figure  50  by  the  receiving  pattern  shown  in  Figure  51.  The 
side  lobe  level  has  gone  down  to  less  than  —50  dB  without  affecting  the  mainbeam. 


Figure  53. 


o  o 

Periodic  DDSA  two-way  pattern  of  signals  from  10  and  -25  related  to 

the  broadside. 
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Figure  54.  Random  subarray  sizes  DDSA  two-way  pattern  of  signals  from  10  and 

o 

-25  related  to  the  broadside. 

D.  SUMMARY  AND  CONCLUSIONS 

In  this  chapter,  the  Three-Step  Matrix  Pencil  method  was  proposed,  and  its 
effectiveness  in  terms  of  RMSE  of  DOA  was  demonstrated.  For  applications  that  require 
a  highly  accurate  DOA  estimation,  the  proposed  method  can  utilize  all  available 
subarrays  to  achieve  high  accuracy. 

Due  to  the  limited  options  for  grating  lobe  suppression  of  a  DDSA  when 
transmitting,  a  combination  of  conventional  methods  plays  an  important  role  in  solving 
this  problem.  By  using  the  combination  of  conventional  methods  on  when  transmitting 
and  the  virtual  filling  method  when  receiving,  an  improved  two-way  pattern  with  an 
ultra-low  side-lobe  level  was  demonstrated. 
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VII.  SUMMARY  AND  CONCLUSION 


This  research  has  focused  on  an  important  and  emerging  issue,  which  is  grating 
lobe  suppression  for  DDSAs.  In  a  general  sense,  using  more  radiation  elements  in  an 
array  system  can  potentially  increase  the  array  gain  when  transmitting  and  SNR  when 
receiving.  However,  because  of  the  mechanical,  structural  and  operational  limitations 
discussed  in  Chapter  I,  the  use  of  a  large  contiguous  array  on  a  platform  can  be  restricted. 
The  idea  of  using  separated  arrays  that  together  form  a  DDSA  is  a  potential  solution  to 
this  dilemma.  The  critical  issues  that  must  be  addressed  are  calibration,  time/frequency 
synchronization,  error  correction  and  grating  lobes. 

This  study  focused  on  various  methods  to  suppress  the  grating  lobes  of  a  DDSA 
on  both  the  transmitting  and  receiving  sides.  Conventional  techniques  used  to  treat  the 
grating  lobe  problem  of  a  single  sparse  array  were  applied  to  DDSAs.  Both  periodic  and 
random  arrangements  were  considered  and  general  equations  for  the  pattern  functions 
were  derived. 

Basic  formulas  for  the  DDSA  radiation  pattern  were  presented  in  Chapter  II.  Both 
periodic  and  aperiodic  or  random  arrangements  were  considered  and  general  equations 
were  derived.  Several  methods  of  potential  interest  for  grating  lobe  suppression  including 
subarray  rotation,  multiplicative  beamforming  and  randomness,  were  discussed  in  this 
chapter  as  well. 

Several  conventional  methods  for  grating  lobe  suppression  were  examined  in 
Chapter  IV.  To  further  improve  the  grating  lobe  suppression  results,  a  combination  of  the 
conventional  methods  was  introduced  and  simulated.  As  can  be  seen  from  the  simulation 
results,  each  method  provides  some  additional  grating  lobe  suppression.  By  combining 
the  methods,  grating  lobes  are  further  lowered  by  around  8  dB  (from  -1 1  dB  to  —19  dB 
for  the  example  presented).  However,  disadvantages  such  as  polarization  loss,  hardware 
complexity,  gain  loss  and  limited  suppression  ability  for  large  separations  restrict  the  use 
of  the  methods,  and  therefore  tradeoffs  need  to  be  made  accordingly. 
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Filling  gaps  between  arrays  with  virtual  elements  for  the  purpose  of  receiving 
processing  allows  a  synthesized  antenna  response  that  duplicates  a  contiguous  array  in  a 
number  of  directions  that  is  limited  by  the  total  number  of  elements  used  in  the 
processing.  The  synthesized  response  has  no  grating  lobes  and  can  have  low  side  lobes, 
as  seen  in  Chapter  V.  As  a  first  step,  the  signal  amplitudes,  phases  and  DOAs  must  be 
extracted  from  the  element  /  and  Q  samples.  This  data  is  used  generate  /  and  Q  samples 
that  would  be  provided  by  virtual  elements  filling  the  gaps  between  distributed  subarrays. 
Low  side  lobes  and  interference  rejection  were  demonstrated  for  the  virtual  processed 
DDSA. 

The  MP  method  was  found  be  to  be  well  suited  for  this  application.  The  MP 
technique  was  extended  to  handle  multiple  subarrays,  for  either  single  or  multiple 
snapshots.  Multiple  snapshots  provide  improved  stability  in  a  low  SNR  situation.  This 
method  can  deal  with  both  coherent  and  non-coherent  signals  and  requires  fewer 
snapshots  for  accurate  DOA  estimation,  as  was  seen  in  Chapter  V. 

The  proposed  filling  method  allows  suppression  of  undesired  signals  that  would 
normally  occur  at  grating  lobe  angles.  It  was  shown  that  for  high  SNR  the  array  response 
approaches  that  of  a  contiguous  array  of  the  same  extent. 

The  Three-Step  Matrix  Pencil  method  was  proposed  in  Chapter  VI,  and  its 
effectiveness  in  terms  of  RMSE  of  DOA  was  demonstrated.  For  applications  that  require 
a  high  accuracy  of  DOA  estimation,  the  proposed  method  can  utilize  all  available 
subarrays  to  achieve  the  highest  accuracy  possible. 

Due  to  the  limited  options  for  grating  lobe  suppression  of  a  DDSA  on  the 
transmitting  side,  a  combination  of  conventional  methods  plays  an  important  role  in 
solving  this  problem.  By  using  the  combination  of  conventional  methods  on  the 
transmitting  side  and  the  virtual  filling  method  on  the  receiving  side,  an  improved  two- 
way  pattern  with  an  ultra-low  side-lobe  level  was  demonstrated. 
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A. 


SIGNIFICANT  CONTRIBUTIONS 


The  primary  contributions  of  this  dissertation  come  from  the  investigation  of  the 
grating  lobe  problem  of  DDSAs  and  the  novel  techniques  proposed  to  deal  with  the 
grating  lobe  suppression  on  the  receiving  side  along  with  the  combination  of 
conventional  methods  to  improve  the  grating  lobe  level  on  the  transmitting  side. 

One  contribution  of  this  research  is  a  proposed  hybrid  approach  that  uses  a 
combination  of  suppression  techniques  on  both  the  transmitting  and  receiving  sides.  The 
result  is  improved  two-way  pattern  performance. 

Another  contribution  is  the  development  of  the  new  virtual  filling  processing 
method  on  the  receiving  side  to  suppress  grating  lobes  and  improve  the  signal-to-clutter 
ratio  and  signal-to-interference  ratio. 

A  final  contribution  is  the  illustration  of  the  relationship  (and,  hence,  tradeoffs) 
between  thermal  noise,  array  errors,  and  the  grating  lobe  suppression  effectiveness.  The 
consideration  of  array  errors  addresses  the  issue  of  array  calibration  and  synchronization, 
which  are  critical  concerns  when  multiple  arrays  operate  coherently. 

B.  FUTURE  RESEARCH 

Future  research  should  focus  on  subarrays  with  different  array  lattices  and  look  at 
the  potential  application  of  the  filling  method.  It  is  relatively  straightforward  to  fill  the 
gaps  for  DDSA  with  the  same  subarray  types  (i.e.,  the  rectangular  periodic  planar  array). 
For  subarrays  with  different  arrangements,  optimal  filling  positions  need  to  be 
determined  for  best  grating  lobe  suppression  performance.  A  similar  question  arises  for 
DDSAs  with  subarray  rotation,  as  well. 

The  mutual  coupling  effect  between  real  array  elements  should  be  included  in  the 
simulation  and  possible  compensation  methods  need  to  be  developed  in  the  future 
research. 

Narrow  band  (single  frequency)  is  assumed  in  this  research.  It  would  be  of  great 
interest  to  extend  the  methods  developed  here  to  broadband  applications.  On  the 
transmitting  side,  broadband  applications  often  assume  small  subarray  sizes.  For  the  same 


87 


number  of  elements  in  a  DDSA,  small  subarray  sizes  will  provide  more  subarrays  or 
subarray  centers;  therefore,  more  randomization  achieved.  On  the  receiving  side, 
parameter  estimation  techniques  for  broadband  applications  need  to  be  examined. 

For  DDSA  design,  parameters  that  will  affect  the  array  performance  are  also  the 
parameters  that  contribute  to  the  grating  lobe  suppression  effectiveness.  For  example,  the 
sizes  of  subarray  gaps,  number  of  elements  in  each  subarray,  element  spacing  in  each 
subarray,  subarray  rotation  angles,  array  element  types  and  array  lattice  are  among  these 
parameters.  Therefore,  the  development  of  an  optimization  algorithm  for  choosing  the 
optimal  parameters  for  best  grating  lobe  suppression  results  will  benefit  the  design 
process  going  forward. 
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